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Neutrino scattering measurements offer a unique tool to probe tiie electroweak and strong inter- 
actions as described by tlie Standard Model (SM). Electroweak measurements are accessible through 
the comparison of neutrino neutral- and charged-current scattering. These measurements are com- 
plimentary to other electroweak measurements due to differences in the radiative corrections both 
within and outside the SM. Neutrino scattering measurements also provide a precise method for 
measuring the 'P2(x,Q^) and x¥z(x,Q^) structure functions. The predicted evolution can be 
used to test perturbative Quantum Chromodynamics as well as to measure the strong coupling 
constant, , and the valence, sea, and gluon parton distributions. In addition, neutrino charm 
production, which can be determined from the observed dimuon events, allows the strange-quark 
sea to be investigated along with measurements of the CKM matrix element \Vcd\ and the charm 
quark mass. 
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I. INTRODUCTION 



Neutrino scattering experiments have strongly influenced the development of particle physics over the past three 
decades. Neutrinos only couple to other particles through the weak interaction, which is well determined in the 
Standard Model (SM). For this reason, neutrino scattering measurements can be used as a probe to measure many of 
the SM parameters or to look beyond the SM for indications of new physics. The effects of the small interaction cross 
section for neutrinos has been overcome by modern experiments through the use of high intensity beams coupled with 
massive detectors which give luminosities in the range of 10"^® cm^^s^^. Data samples in excess of one million events 
are now available which allow measurements of strong and electroweak parameters comparable in precision to other 
fixed-target and collider determinations. 

This Review addresses three broad physics topics accessible to high energy neutrino d eep inelastic scattering exper- 
iments: nucleon structure functions and tests of quantum chromodynamics (QCD) (Sec. |lll| ), neutrino induced charm 
production (Sec. iV), and electroweak physics with neutrino beams (Sec. M). We begin with a brief introduction to 
the relevant kinematics and cross section formulae, and follow with a summary of experiments whose results we review 
(Sec. A| ). We then devote an entire section to the difficult problem of normalizing data through the determination 
of the neutrino flux and total cross section (Sec. |l|). The main physics sections follow. 

Experiments we review date from roughly 1980 until the present, although in some cases we reach back further if 
the results are still of interest. We also restrict ourselves to high energy measurements; thus we emphasize the CERN 
SPS and Fermilab Main Ring and Tevatron programs. Earlier lower energy experiments have been reviewed elsewhere 
(Barish, 1978; Fisk and SciuUi, 1982; Diemoz et ai, 1986; Mishra and SciuUi, 1989) 



A. Kinematic Formalism 




FIG. 1. The first order Feynman diagram for deep inelastic neutrino scattering 



The tree-level diagram for charged-current (CC) neutrino-nucleon (i^N) scattering is shown in Fig. |l|. A neutrino 
(antineutrino) with incoming four-momentum ki scatters from a quark or antiquark in the nucleon via exchange of 
a (^ ) boson with four-momentum q. In the laboratory frame, the variables that can be measured in this 
interaction are the energy and angles of the outgoing muon, E^, 6^, (f>^, and the energy and angles of the outgoing 
hadrons Ehad, Shad, (f'had- In practice, experiments are located sufficiently far from the neutrino production point that 
the neutrino beam can be considered, to high accuracy, parallel to the z-axis, thus the difhcult-to-measure hadron 
scattering angles are not necessary to specify the kinematics. In the laboratory frame, the 4-vectors of the event can 



then be written in terms of E"^, 0^ , 0^, Ehad, and M, the proton mass:a 



fci = {E,,0,0,E,), 

h = {Ef,, Ef,sin9^cos4)^, Ef,sm9f,siii(f>^, Ef,cos9^), 
p=(M, 0,0,0), 
q = ki- k2, 



(1) 



^We use the conventional units of high energy physics in this Review: h — c — 1, and all dimensionful quantities are in GeV 
unless otherwise indicated. 
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with E,, = E„ 



Eh 



Useful invariant quantities which describe the interaction are: 



In the laboratory, these reduce to: 



^ = (p ■ q)/M (energy transfer), 
y — P ' q/p ■ ki (inelasticity), 

= —q^ (negative squared 4 — momentum), 
X ^ / {2p ■ q) (Bjorken scaling variable). 

= Ehad - M, 
[Ehad - M) 



(2) 



{Ehad - 
{Ehad - 
{Ehad 



E,i)Ef_,9f^ 



(3) 



2M{Ehad-M) 

using the small-angle approximation and neglecting the muon mass. 

In principle, for a given beam of neutrino energy E^, the accessible kinematic region is bounded by: 



m. 



%xMEl 

4 

^xMEl 

4 

mf 



<y< 



< 



E. 


(1 


+ 2Mx/E^) 
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(1 


+ 2Mx/E^) 




2ME^x 


(1 


+ 2Mx/E^) 



E, 



(4) 



2ME„x. 



<x<l. 
2ME^ - - 

The lower limits are given to leading order in the final state lepton mass and can always be set to for the high 
energies considered; the arrows indicate the high energy limits of the upper range on the variables, which, for most 
purposes, are also sufficient. In practice, experimental restrictions to select regions of high acceptance and small 
corrections reduce the available range. The accessible kinematic range of a given experiment can be described in the 
plane of any two of the above variables. The resulting kinematic range in the x and plane is shown in Fig. ^ for 
four high-statistics deep inelastic neutrino scattering experiments. 
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FIG. 2. The kinematic regions accessible to the high statistics neutrino experiments discussed in this Review. 



In the case of neutral current (NC) scattering, the scattered neutrino is not reconstructed; all event information 
must be inferred from the hadron shower. The important kinematic quantities are identical to the charged-current 
case up to small effects due to the muon mass. 



B. The Neutrino Scattering Cross Section 

The unpolarized cross section for neutrino scattering is derived from placing the constraints of Lorentz invariance 
and CP invariance on the product of the leptonic and hadronic tensors. The result is that, for fixed x and Q^, only 
five free parameters describing the hadronic structure remain. Two are proportional to the outgoing lepton mass, 
which is small for the case of scattering and can be dropped. The remaining three parameters are functions of x 
and and are called "structure functions." 

In this Review, we will describe the cross section for scattering from an isoscalar target in terms of the structure 
functions F2, xF^ and Rl defined by: 
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dxdy TT (1 + Q^/M^) 

■ . / + (2Mxy/Qr ^,^y_Mxy\ 

±.F3^(-)%(l-|)]. (5) 

where the ± is +(— ) for 1^(17) scattering. Analogous functions F2 '^■'^(a;, Q^) and Q^) appear in the cross 

section for deep inelastic charged lepton scattering. Rl{x, Q^) can be interpreted as the ratio of the longitudinal to 
transverse virtual boson absorption cross-section, and is related to the 2x^1(0;, Q^) structure function by: 

^ F2{x, Q^){l+4APxyQ^)-2xFi{x, Q^) 

2xFi{x,Q^) ■ ^ ' 

The function xFz{x,Q^) appears only in the cross section for the weak interaction. It originates from the parity- 
violating term in the product of the leptonic and hadronic tensors. 

C. Neutrino Scattering from Partons 

The cross section presented in Eq. ^ makes no assumptions about the underlying structure of the hadron involved in 
the interaction. However, the results of early experiments provided clear evidence that the neutrinos were interacting 
with quarks in the nucleon (Perkins, 1975). This section discusses neutrino-parton scattering at the simplest level. 



The interpretation of the cross section within QCD is discussed in more detail in Sec. [II. 

The "naive" quark-parton model (QPM) describes the structure functions in terms of momentum-weighted parton 
distribution functions (PDFs). For neutrino interactions with the partons in the proton, 

F2!hpMi^) = 2a: id{x) + u{x) + six) + c{x)) , (7) 
■^^a'QPAfC^) ~ 2x {d(x) — u{x) + s{x) — c{x)) . (8) 

In these equations u, d, s and c indicate the PDFs for up, down, strange and charm quarks, respectively, in the proton. 
These PDFs describe the probability that the interacting parton carries a fraction, x, of the proton's 4-momentum.El 
Note that because the neutrino interacts via emission of a W~^, only the negatively charged quarks contribute to the 
interaction. 



^This interpretation of x is an approximation that we employ due to its intuitive appeal, and because at high energies, this 
definition of x differs negligibly from its more rigorous interpretation as a light-cone- variable (See Sec. [[ A|.). 
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In a similar way, neutron structure functions can be written in terms of the proton PDFs by invoking isospin 
invariance, 

F^Qpnix) = 2a: {u{x) + d{x) + s{x) + c{x)) , (9) 

^P3.qpm{x) = 2x {u{x) - d{x) + s{x) ~ c{x)) . (10) 

When neutron and proton structure functions are combined for an isoscalar target, one obtains: 

F!^Zpm{x) = X iH^) + d{x) + 2s{x) (11) 
+u{x) + d{x) + 2c{x)) , 

"3,QPAf ^ 

—u{x) — d{x) — 2c(a;)) 



xFjqpm{x) = X {u{x) + d{x) + 2s{x) (12) 



Antineutrino interactions proceed via emission of a W~ boson. Therefore only the positively charged quarks 
contribute. Following the prescription for the neutrino case and assuming that s{x) = s{x) and c{x) = c{x), the 
antineutrino structure functions are: 

F2,QPMi^) = F2^Qp]^f{x) (13) 

xFl^PMix) - xFl'^p,,ix) - Ax (six) - c{x)) . (14) 

In this review xF^, = ^{xF!^'^ + xF!^'^) will refer to the average of the neutrino and antineutrino measurements; 
^{xF^) = xF^^ — xF^^ will indicate the difference between the neutrino and antineutrino structure functions. We 
use F2 to refer to either neutrino or antineutrino scattering because the structure functions are equivalent. 

In the QPM, the interaction probability for each PDF depends on the square of the charge associated with the 
interaction. In the weak interaction, this is unity. For charged-lepton scattering mediated by a virtual photon: 

^2,qpm(2;) = i^i^) + d{x) + u{x) + d{x) 

+ '^{s{x) + -s{x)) + ^{c{x)+c{x)) . (15) 

The structure function definitions of the QPM must be modified to accommodate strong interactions between the 



partons and to include mass effects. Within the framework of QCD, discussed in Sec. [II A 1, the parton distributions 
acquire a dependence. Small transverse momentum quark-giuon interactions may occur at the time of the scatter. 
The probability that the interaction resolves such processes depends on and modifies the parton distributions, 
leading to a dependence known as "scaling violations." Thus at Leading Order (LO) in QCD, the structure 
functions are: 

F2,LO= J2 xqix,Q^)+xq{x,Q^), (16) 



xF 



3,LO 



E 



xq{x, Q^) - xq{x, Q^). 



(17) 



(In the above equations and subsequent sections, q and q are introduced as the generic notation for the various 
parton distributions in the proton.) Parton mass effects further modify the parton distributions, an effect called "slow 
rescaling" (see Sec. iVB), which leads to q{x, Q^) q{£,,Q^), where ^ depends on the mass of the parton. 

To LO, RLix,Q'^) = and the Callan-Gross relation (Callan and Gross, 1969), F2{x,Q'^) ~ 2xFi{x,Q'^), is valid 
because spin-i massless quarks cannot absorb a longitudinally polarized virtual boson, reverse direction in the center 
of momentum frame, and still conserve helicity. A non-zero Rl{x,Q'^) is a manifestation of non-perturbative and 
higher order QCD effects and is described in Sec. [II A 1. In all formulations that we are aware of, the structure 



function i?^ is the same for neutrino and antineutrino scattering. 

These tree-level expressions can be used to point out some of the interesting information accessible from vN 
scattering: 



F2 (x, Q ) and F2 [x, Q ) measure the sum of quark and anti-quark PDFs in the nucleon. 
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The average of xF^^ {x, Q^) and xF^^{x, Q^) measure the valence quark PDFs, 



(18) 



F2^ {x,Q^) and F!^^'^^'^ {x,Q'^) together confirm the fractional electric charge assignment to the quarks and 
provide sensitivity to the strange quark PDF. 

xF^^ {x,Q^) — xF^^ {x,Q^) is sensitive to the strange and charm quark content in 



• The difference A(a;F3) 
the nucleon. 



Neutrino scattering from d or s quarks can produce a c quark in the final state which can be identified by either a 
hfetime or, more easily, a semileptonic decay tag. This ability to tag final state charm permits a series of unique strong 
and electroweak tests including the direct measurement of s{x,Q'^) and s{x,Q'^), extraction of the CKM elements 
Vcd and Vcs, a measurement of the mass of the charm quark rric, and a sensitive probe of the dynamics of LO and 



next-to-leading order (NLO) QCD. This physics is described in Sec. |V. 

Neutral-current vN scattering is described by a closely related set of structure functions. While these can be 
measured, the primary interest in this process lies in inclusive cross-section measurements for electroweak physics 
tests of the Standard Model. The detailed description of i^N NC scattering will be discussed in Sec. 



D. Overview of Experiments 

This section provides an overview of the beams and experimental apparatus used in deep inelastic neutrino scatter- 
ing. Table || summarizes the data samples obtained by the high-precision neutrino experiments which are discussed 
in more detail below: 



Experiment 


Beam Energy (GeV) 


target 


ly events 


1/ events 


Flux Method 


CCFR (E616) 


30-300 


Fe 


150,000 


23,000 


Secondary Monitoring 


CCFR (E701) 


30-300 


Fe 


35,000 


7,000 


Sec. Mon., y-int. 


CCFR (E744/770) 


30-360 


Fe 


1,300,000 


270,000 


Fixed 1/0 


CDHSW 


20-212 


Fe 


221,000 


20,500 


iter, y fit 


CHARM 


10-160 


CaC02 


116,000 


6000 


Secondary Monitoring 


CHARM H 


5-100 


Glass 


750,000 Q.E. 


1,300,000 Q.E. 


"tot 


BEBC 


15-160 


Ne (D2) 


15,000 (12,000) 


10,000 (11,000) 


iter, y fit 



ref 



(MacFarlane el al, 1984) 
(Auchincloss, et al., 1990) 
(Sefigman, 1997) 
(Berge et al, 1987) 
(Allaby et al, 1988) 
(Geiregat et al., 1993) 
(Allport et al, 1989) 



TABLE I. Overview of High Precision Neutrino Experiments. Target refers to the bulk of the material in the experimental 
target. Total events vary for each type of analysis depending on cuts; the numbers here are approximate. 
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1. CCFR Experiment 



The CCFR collaboration has performed a series of experiments at Fermi National Accelerator Laboratory (hereafter, 
"Fermilab") using the Lab E neutrino detector. This paper includes results from data runs E616, E701, E744 and 
E770, spanning a period from 1979 through 1988. A new experiment, E815, is presently taking data. 

The E616 and E701 experiments took data using the Fermilab dichromatic neutrino beam. The primary protons 
had an energy of 400 GeV. The secondary beam traversed a 60 m chain of magnets and collimators, which selected 
pions and kaons based on energy and charge. The neutrino energy spectrum from this beam has two relatively narrow 
peaks from the pion and kaon decays which gives the beam its dichromatic nature. 

The most recent results from the CCFR CoUaboration come from the 1985 (E744) and 1987 (E770) fixed-target 
runs at Fermilab. The CCFR neutrino beam resulted from decays of pions and kaons produced in interactions of 
the 800 GeV FNAL proton beam with a beryllium target. The FNAL Quadrupole Triplet beam line, which had no 
sign-selecting magnets, was used to transport the secondaries. This resulted in a wide-band beam of neutrinos and 
antineutrinos, with a relatively uniform energy spectrum ranging in energy up to 600 GeV. This beam resulted in 
approximately 86.4% v^^, 11.3% and 2.3% Ve and Ve events. 



Blue Cart 




FIG. 3. Schematic representation of the CCFR detector. The neutrino beam travels from left to right. The target-calorimeter 
is on the left and the muon spectrometer (toroid) is on the right. 

Neutrino events were observed in the Lab E neutrino detector (Sakamoto, et ai, 1990), shown in Fig. ^. The detector 
is constructed as a target calorimeter followed by a toroid muon spectrometer. The 690 ton calorimeter consists of 
eighty-four 3 mx3 mxlO cm iron plates interspersed with scintillators and drift chambers, with a measured hadronic- 
energy resolution a/Ehad — 0-89 /^/Ehad- The toroid spectrometer has five sets of drift chambers for muon tracking as 
well as hodoscopes for triggering; its momentum resolution is limited by multiple Coulomb scattering to op/p = 0.11. 
The detector was periodically calibrated by a test beam of pions and muons. 

NuTeV (Bolton et ai, 1990), the next generation of deep inelastic neutrino scattering experiments at Fermilab, 
began taking data in May, 1996. NuTeV has modified the beam line to select the sign of the charged pions and kaons 
(Bernstein et ai, 1994). The resulting beam is almost purely v^, or F^, depending on the sign selection. The enhanced 
purity of the antineutrino sample will improve a number of QCD studies, especially those concerned with small- a; 
physics. Separate running will also permit the first high-statistics simultaneous measurements of the electroweak 
mixing parameter sin^ 9^ and neutral-to-charged-current coupling strength ratio p. The substantially upgraded Lab 
E neutrino detector is continuously calibrated with a test beam of pions, electrons or muons throughout the running. 
The new test beam also provides a broader range of hadron and muon energies than were available to CCFR. 

2. The CDHSW Experiment 

The CDHSW experiment (Abromowicz et ai, 1981; Holder et ai, 1978) measured the total neutrino cross section 
using 100, 160 and 200 GeV narrow band neutrino beams and performed precision electroweak and structure function 
measurements with a wide-band beam during the early 1980's. Primary protons from the CERN SPS were used to 
produce secondary pions and kaons. For the wide-band running, the secondary beam was focussed using a magnetic 
horn. Depending on the sign selection set by the magnetic field, the tertiary beam consisted mainly of neutrinos 
or antineutrinos. An antineutrino contamination in a neutrino beam will result in single-muon events with a muon 
that has positive rather than negative charge. This contamination is therefore called the "wrong-sign" background. 
The fraction of wrong-sign background events for CDHSW neutrino running was 2.3%. For antineutrino running, the 
wrong-sign background was 17% (Berge et ai, 1991). 

This experiment used the upgraded CDHS detector, which consisted of toroidally-magnetized iron plates sandwiched 
between planes of scintillator and drift chambers. The calorimeter was 3.75 m in diameter, 21 m long, and had a 
total mass of approximately 1100 tons. The hadronic energy resolution of the detector varied, depending on the 
segmentation of the iron, from cr/Ehad ~ 0.58/^/Ehad to 0.70/^/Ef^. The momentum of the muon was determined 
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by the bend in the magnetic field, with an average resolution of 9%. The field of the toroid was chosen to focus 
or fjT based on the setting of the magnetic horn. 



3. The CHARM Experiment 



The CHARM experiment (Jonker et aL, 1982) made a high precision measurement of the total cross section for 
neutrinos and antineutrinos using the 160 GeV narrow band neutrino and antineutrino beam at the CERN SPS 
during the early 1980's. The CHARM detector was composed of a fine-grained calorimeter followed by a muon 
spectrometer. The calorimeter consisted of 78 marble plates (CaCOs) of dimension 3 mx3 mx 8 cm. Scintillation 
counters, proportional drift tubes, and streamer tubes were located between the plates. The resulting hadronic 
energy resolution was a / Ehad = O-'i^/y/Ehad- The vertex position resolution was 3 cm at Ehad = 50 GeV. The muon 
spectrometer consisted of an iron toroid spectrometer surrounded by a magnetized iron frame and instrumented with 
proportional chambers. The calorimeter was calibrated with electron and pion beams. 



4. The CHARM II Experiment 



The CHARM II experiment (Geiregat et al, 1993) collected data during the period 1987-1991, using the wide-band 
neutrino beam described in the above CDHSW experiment section. This experiment was optimized for detection 
of electrons. The detector's target was composed of 48 mm(0.5 radiation length) thick plates of glass, with a total 
mass off 690 tons. Interspersed between the plates were streamer tubes, and between every five sets of plates and 
streamer tubes were scintillation counters. The electromagnetic energy resolution from the streamer tubes was a/E = 
0.09 -1-0. 15/-v^. The muon spectrometer consisted of magnetized iron toroids instrumented with scintillation counters 
and drift chambers. The resolution was AP/P = 13% at 20 GeV. 



5. The BEBC Experiment 



The Big European Bubble Chamber (BEBC) experiment (Baruzzi et al., 1983) ran with both narrow- and wide-band 
beams provided by the CERN SPS during a series of experiments that began in the late 1970's and continued during 
the 1980's using a variety of targets. This Review mainly discusses results from the high statistics BEBC runs with 
deuterium (WA25) and neon (WA59) targets. The cryogenic bubble chamber had a 10 m^ fiducial volume surrounded 
by a 3 T superconducting magnet and supplemented with a two-plane, 150 m^ external muon identification system. 



II. NEUTRINO FLUX AND TOTAL CROSS SECTION MEASUREMENTS 



Accurate measurement of the flux of neutrinos in a beam is necessary for all of the analyses discussed in this Review. 
Flux monitoring is a difficult task for a beam of neutral weakly interacting particles; hence indirect methods must 
be employed. Flux determinations can be based on measurements of the secondary beam. Alternatively, experiments 
have used fits to the y distribution of the data or integration over the observed low y events. These methods rely 
upon the relation between the cross section, flux and number of events (A^i/,f): 



and on assumption that as y — > 0, 



1 r/rr"^ 1 da"^ 



Ev dy El, dy 

where C is a constant independent of energy. The Flux measurement errors are among the largest systematic errors 
associated with the neutrino precision measurements, therefore an extended discussion of the subject is presented 
here. 
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A. Absolute Flux Determination 



Neutrino beams result mainly from the decay of secondary pions and kaons which are produced when a proton 
interacts with a target. Most pion and kaon decays are two-body, thus the neutrino energy and flux can be determined 
using kinematic constraints if the integrated intensity, energy, and spacial distribution of the secondary beam are 
known. This method is most successful when applied to narrow-band beams and was employed by many experiments 
to measure structure functions and the total cross section. These experiments include CCFR E616 (Oltman et al., 
1992; Auchincloss, et at, 1990), CDHSW (Berge et at, 1987), and CHARM(Allaby et ai, 1988). 

The primary and secondary beams are monitored using a variety of detectors, including RF cavities, toroids, 
Cerenkov counters, segmented wire ionization chambers, and hodoscopes. Redundancy is important to reduce the 
systematic errors in the measurement of the secondary beams. Two uncertainties associated with this technique 
result from three-body decays of kaons and pion or kaon decays that occur before the momentum selection magnets 
(Auchincloss, et ai, 1990). The former are largely from K ir^i', which can be corrected using a Monte Carlo 
simulation of the neutrino flux. The latter correction can be measured directly by placing a dump in the beam 
after the magnet, thus absorbing the sign-selected secondary beam. All interactions in this mode would result from 
neutrinos produced in decays before the magnet. As an additional constraint to reduce errors further, the CHARM 
collaboration also monitored the muon flux from the tertiary beam (Allaby et al., 1988). 



B. Relative Flux 



The structure function measurements described in Sec. Ill are very sensitive to the accurate determination of the 
relative flux among the various energy bins and between neutrinos and antineutrinos. Therefore, it is important 
to establish methods of normalizing between neutrino energy bins which have small systematic errors and little 
dependence on assumptions of nucleon structure. These methods are crucial for experiments using wide-band neutrino- 
beams, where beam monitoring methods cannot provide sufficient accuracy. In a narrow-band beam, these methods 
can be used to test the accuracy of the measurements from secondary beam monitoring. 

The "y-intercept technique" relies on the premise that near zero hadronic energy transfer, y = Ehad/Ei, — s- 0, 
the differential inelastic cross section divided by energy is independent of energy and is the same for neutrinos and 
antineutrinos (Eq. |20|).The differential cross section for neutrino scattering, given by Eq. |^, in the y limit reduces 
to: 



- 1 


da"- 




■ 1 


da"- 






dy 


y=0 




dy . 


y=o 



^— / F^ix.Q^ ^Q)dx^ Const. 



(21) 



This is true independent of incident neutrino energy or neutrino type for scattering on an isoscalar target. Thus, for 
each bin in energy, the flux of neutrinos or antineutrinos is given by: 



lim , 

y^o V Ev 



1 dW^''^{E^) 



dy 



Const. X (£;^). 



(22) 



This method was applied by CCFR (E616) to cross-check the flux obtained from secondary-beam measurements for 
the determination of the neutrino total cross section. 

For wide-band neutrino beams, the beam energy and radial distribution arc usually determined through iterating 
the Monte Carlo simulation of the neutrino events, successively varying the input flux and structure functions until 
the Monte Carlo matches the observed distribution of events. Initial distributions for the structure functions are taken 
from previously published results. These structure functions and initial distributions for the neutrino flux are used 
to generate events which are then put through the Monte Carlo simulation of the detector. The ratio of simulated 
events to detected events in each x-y-E^, bin provide the correction factors for the flux and structure functions for the 
next iteration. Variations on this method are used by both the CDHSW and CCFR collaborations. 

The CCFR E744/E770 "flxed-z^o method" provides a specific example of one iterative technique for extracting the 
relative flux (Seligman, 1997). This method exploits the y-dependence of the data to obtain the flux. The relationships 
of Eq.'s 19 and can be extended to small as a polynomial in terms of ?/ = v / E^^ 

dN^,i7 



dv 



= A,j, + B,j,{v/E,) + {C,..vl2){v/E,Y, 



(23) 
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where A, B and C are given by 

A 
B 



GlM 



F2{x,Q^)<^{E,)dx, 



Jo 



(24) 



C = B- 



GlM 



with 



F2{x,Q^)R{x,Q^)^{E,)dx, 



1 + 2Mx/v Mx 



l+i?L(.T,Q2) y 



(25) 



Except for small variations due to the scaling violations of the structure functions, A, B and C are relatively 
independent of u and Ei,. Integrating equation 23, the number of events with v < vi^ for a given energy, N^^1^{Ei,)^ 
is given by: 



N^XiE,)^m.) dvA 



V B v"^ ( B f F2R 
E^A~2E^, ll^TTT 



(26) 



From the values of N^^l^^{E^), $(£;^) can be determined for both neutrinos and antineutrinos up to an overall 
norm alization constant. The overall normalization is determined from to the world-average cross-section for vFe (Sec. 
ilD|) (Seligman, 1997), 

a^^yE^ = (0.677 ± 0.014) x lO^^^cmVCeV. 

The assumption that the integrals over the structure functions are independent of v is only approximately true. 
Recall that = 2Mvx, therefore for fixed using F2 as an example, 



(27) 



F2{x, Q') = / F2{x, 2Mvx) = F^^\v) 



The systematic error introduced by this residual i/-dependence was studied by varying the v range from which the 
flux was extracted and found to be less than 0.5% (Seligman, 1997). 



C. Electron Neutrino Flux 



Most neutrino beams are composed mainly of muon-type neutrinos and antineutrinos, with a small contamination 
from the electron flavor. Accurate determination of this contamination is important for the precision measurement of 
sin^ 9w discussed later. The Vf./ve flux is evaluated through detailed Monte Carlo simulations which are tied to the 
direct muon- flavor flux measurements. As an example of the magnitude of the contamination, the CCFR E744/770 
Monte Carlo Flux determination is shown in Fig. ^ (King, 1994). As a test of the accuracy of this method, a recent 
analysis of the CCFR data to determine directly the total number of Ve interactions has confirmed the Monte Carlo 
prediction (Romosan et ai, 1996). 
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FIG. 4. The CCFR 744/770 neutrino and antineutrino event spectrum. 



D. Measurements of the Total Cross Section 



The total neutrino and antineutrino cross sections have been measured in a large number of experiments. This 
section reports on those experiments which have measured the total cross section to an accuracy of at least 10%. 
The CHARM (Allaby et ai, 1988), CDHSW (Berge et ai, 1987), BEBC (Aderholz, et al., 1986; Bosetti et at, 1982), 
CCFR E616 (MacFarlane et ai, 1984), and CCFR E701 (Auchincloss, et ai, 1990) results are from experiments 
described in detail above. Two older experiments, the 15-ft Bubble Chamber (15-ft BC) at Fermilab (Baker et ai, 
1983; Taylor et ai, 1983), and FNAL E310 (WHPFOR), a counter experiment with an iron target (Heagy et al, 
1981), also have results with small errors and are discussed below. All results are presented for interactions with an 
isoscalar target, and a correction has been applied to the iron data to compensate for the neutron excess. 

A linear dependence of the total cross section on energy is expected if point-like scattering between neutrinos 
and quarks dominates the scattering mechanism. This behavior is observed in the neutrino-iron total cross section 
measurements from CCFR E616 (MacFarlane et al, 1984), CCFR E701 (Auchincloss, et ai, 1990) and CDHSW 
(Berge et ai, 1987), as shown in Fig. |^. The ratio (Jtot/Ey as a function of E,^ is in agreement with a straight line 
with no slope for both neutrino and antineutrino interactions. 



Ll/Fe 
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FIG. 5. Neutrino (solid) and antineutrino (open) measurements of atot/E^. Left plot: atot/Ev versus E^ for iron targets. 
Circles - CCFR E616, Squares- CCFR E701, triangles - CDHSW. Right plot: average atot/E for a variety of targets (isoscalar 
corrected). Dashed lines indicate average for iron data (see text). (The error bars include both statistical and systematic 
errors.) 
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The total cross section for neutrino scattering on iron is used to provide the overall normalization of the flux 
for the CCFR E744/E770 experiment. Table summarizes the isoscalar corrected measurements for neutrino-iron 
experiments. The world averages for these experiments are indicated by the dashed lines in Fig. ||(a). The data 
from all experiments arc in good agreement with these averages. The E744/E770 data also yield a measurement of 
a~/a''=0.509 ± 0.010 (Scligman, 1997) which is listed with the other measurements in Table ||. 



Experiment 




a~/E 


a~/a" 


ref. 


CCFR (E616) 


0.669 ± 0.024 


0.340 ± 0.020 


0.499 ± 0.025 


(MacFarlane et ai, 1984) 


CDHSW 


0.686 ± 0.020 


0.339 ±0.022 


0.495 ±0.010 


(Berge et ai, 1987) 


CCFR (E701) 


0.659 ± 0.039 


0.307 ±0.020 


0.467 ± 0.028 


(Auchincloss, et al., 1990) 


CCFR (E744/770) 






0.509 ±0.010 


(Seligman, 1997) 


World Ave. 


0.677 ±0.014 


0.334 ±0.008 


0.500 ± 0.007 





TABLE IL Measurements of the neutrino cross section on iron 
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FIG. 6. Giy/ Ey for and from CCFR E744/E770. Normalization is based on the world averages given in Table ^ 



The world-average v/v cross-section ratio, in combination with the world-average neutrino cross section, sets the 
overall normalization for CCFR E744/E770 data. Given this normalization, CCFR has measured the energy depen- 
dence of a tot I up to energies of 400 GeV. The result is shown in Fig. where systematic and statistical errors 
have been added in quadrature. The data have been fit to a model which allows a slope in UtotjEi,^ denoted ei,(p) for 
neutrinos (antineutrinos), 

arifV£;. = a'^(^)(l + (28) 

the results are 

= (-2.2 ± 0.7)%/(100 GeV) (CCFR), 
£p = (-0.2 ± 1.3)%/(1GG GeV) (CCFR). 

A small slope for the neutrino cross section is expected as a consequence of QCD and heavy quark effects (Hinchcliffe 
and Smithy 1977). 

Figure |5|(b) shows the average total cross section measurement for neutrino-iron targets as well as neutrino-nucleon 
scattering from various isoscalar targets by CHARM (AUaby tt ai, 1988), the 15-ft BC (Baker et ai, 1983; Taylor 
et ai, 1983), and BEBC (Aderholz, et ai, 1986; Bosetti et al., 1982). The approximate average beam energy of 
these lower energy measurements is 25 GeV. The results are in good agreement with the average for the high energy 
neutrino-iron cross section indicated by the dashed line. 



III. MEASUREMENTS OF NUCLEON STRUCTURE FROM INCLUSIVE SCATTERING 

Precision deep inelastic neutrino scattering experiments provide opportunities to test QCD evolution and extract 
the QCD parameter A, which sets the scale of the strong interaction. In the kinematic regions where the structure of 
the nucleon can be interpreted in terms of quarks, neutrino scattering possesses high sensitivity to many individual 
parton distributions. Probing the nucleon with neutrinos complements charged-lepton scattering experiments, and 
comparisons between neutrino and charged-lepton deep inelastic experiments provide tests of the universality of the 
structure functions. Global analyses which include the neutrino data and other cross section measurements yield the 
best parameterizations of the individual parton distributions. 

This section reviews the contributions of neutrino deep inelastic scattering experiments to precision measurements 
of nucleon structure. A brief review of the theory of the structure functions is provided. Experimental techniques and 
issues which must be considered in the analyses are reviewed. Results of the high precision neutrino measurements are 
compared between experiments and to charged-lepton scattering results. Finally, the QCD analysis of the structure 
functions is presented. 



A. Theoretical Framework 



This section describes the theoretical framework in which the structure function measurements will be interpreted. 
The QCD expectation is described along with a discussion of non-perturbative and nuclear effects. 



1. The Structure Functions within QCD 

The general form for the differential cross section depends upon three structure functions, F2, Rl and xF^^ as 
shown in Eq. |[ In the QPM, the neutrino structure functions can be written as functions of sums and differences 
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of the momentum- weighted quark probabihty densities, as described in Sec. iC. QCD modifies this parton model 
interpretation in order to account for the interactions between the partons. 

Within QCD, partons with higher fractional momentum may contribute to interactions at any lower x through 
radiation or gluon splitting. The probability that the interaction resolves such a splitting depends on Q^, as described 
by the Dokshitser-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations (Altarelli and Parisi, 1977; Dokshitser, et ai, 
1978; Dokshitser, et ai, 1980; Gribov and Lipatov, 1972): 







dlnQ2 


27r 


dq^{x,Q^) 




din 02 


2tt 




+ G{y,Q' 


dG{x,Q^) 


asm 


dlng2 


27r 



(29) 



y 

^ f^-^[q'{v,Q')P^,{x/y) 

^ J X y 

i,Q^)P,G{x/y)] , 

Jx y 

+ G{y,Q^)PGGix/y)] . 

The splitting functions, Pij{x/y), give the probability that parton j with momentum y will be resolved as parton i 
with momentum x < y . The symbols q^^ = "^{qi — qi) and q^ = '^{qi + qi) refer to the non-singlet and singlet 

i i 

quark distributions, respectively. The probability of finding a gluon in the nucleon carrying a fractional momentum 
X is represented by G{x, Q^). 

From the DGLAP equations, it can be seen that the change in the structure functions with as a function of x 
depends on a^. At next to leading order, as is given by: 



47r 



/3oln(g7A2) 



1 



A ln[ln(QVA2)] 
(3o In(g2/A2) 



(30) 



where (3o — 11 — 2nf /5 , (3i — 102 — 38n//3, and nf is the number of quark flavors participating in the interaction at 
this Q^. The QCD coupling constant A must be introduced when the renormalization technique is applied to remove 
the divergences within QCD. For the discussion here, the modified minimal subtraction scheme {AIS) is used. 

The structure function Rl = (yL/o'T is zero in the simple parton model. However, QCD effects such as quark-gluon 
bremsstrahlung and qq pair- production introduce transverse momentum leading to a small, non-zero value oi Rl 
which has a predicted x and dependence. Within QCD (Altarelli and Marinelli, 1978): 



Rl,qcd{x,Q'^) 



(31) 



IT J X 



1 dy 



E q.ix)a^ + J: Gix)a 



%=q,q 



E q.{x)a^ (f ) + E G{x)a^ U 

i=q,q i=q 



Given well- measured quark distributions and as from xFj, and F2, a precise measurement of Rl can be a sensitive 
probe of the gluon distribution. 



2. Nonperturbative QCD Effects 

Non-perturbative QCD processes which contribute to the structure-function measurements are collectively termed 
higher-twist effects. These effects occur at small where the impulse approximation of scattering from massless non- 
interacting quarks is no longer valid. Examples include target mass effects, diquark scattering, and other multiparton 
effects. Because neutrino experiments use heavy targets in order to obtain high interaction rates, nuclear effects must 
also be considered. 

a. Target-Mass and Higher- Twist Effects. The target-mass correction (Georgi and Politzer, 1976) to the structure 

functions accounts for the mass of the nucleon M by rescaling the apparent fractional momentum of the quark: 

1 

X ^ ^ = 2x1(1 -t- fc), where = (l -I- 'Ix^M'^ jQ'^^ ^ . The structure functions are then: 
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TMC 



In the above equations, 



^pTMC 



h 
h 
h 



QCD 



fc2 ^2 
fc2 ^2 



6M^ x-^ 

"Q2~fcl 1 

x^ 



(32) 



dvF: 



QCD 



(«,Q2)/„2 



(33) 



= duxF^^^{u,Q^)/i 



For the QCD analysis presented in Sec. HIE, the structure functions are corrected to remove the target- mass effect. 

The remaining higher-twist effects cannot be calculated a priori from PQCD, but must be measured. Virchaux 
and Milsztajn (1992) compared SLAC electron scattering and BCDMS muon scattering measurements of F2 to QCD 
expectations. Deviations were attributed entirely to higher-twist effects and were fit to the form: 



F 



Measured 



edicted 



(1 + C,/Q^). 



(34) 



The constants, C^, calculated for each a;-bin are shown in Fig. 0. This determination was done using a value of 
smaller than the current world average which may make these higher twist corrections too large. 
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FIG. 7. Coefficients of the Higher-Twist parameterization. Solid - xFs from Sidorov (1996). Open - F2 from Virchaux and 
Milsztajn (1992). The solid(dashed) line is the F2{xF3) result from Dasgupta and Webber (1996). 

The Virchaux and Milsztajn (1992) corrections are applicable to a measurement of F2 in neutrino scattering, 
however xFt, may not have the same higher-twist corrections as F2 . Figure ^ also shows the measured higher- twist 
contribution extracted from the CCFR xF^ data based on a NLO analysis following the prescription used for F2 
(Sidorov, 1996). In addition. Fig. ^ shows the calculation by Dasgupta and Webber (1996) which includes infrared 
renormalon contributions. The measured higher-twist contribution to the two structure functions are consistent; 
however, the coefficient for the x = 0.65 bin of xF^ is almost a factor of two lower than that for F2. A preliminary 
analysis of the recent CCFR F2 and xF^ results also indicates higher-twist corrections that are smaller than the 
Virchaux and Milsztajn (1992) or Dasgupta and Webber (1996) determinations. For these reasons, the new CCFR 
analysis uses one-half of the Dasgupta and Webber (1996) values, with a systematic error given by repeating the 
analysis with no correction and with the full correction. 

b. Nuclear Effects in Neutrino Scattering Dependence on the mass number A of the nucleus could arise from 
several effects: 

• Fluctuations of the intermediate virtual boson to mesons. This might suppress the bound nuclcon (iron) structure 
function compared to that from a free (hydrogen) or nearly free (deuterium) nuclcon (Bauer et al, 1978). The 
vcctor-meson-dominance model ascribes the cause of shadowing to fluctuations of the boson into mesons leading 
to strong interactions near the "surface" of the nucleus. This model was developed for charged-lepton scattering. 
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where the photon can fluctuate only into vector mesons. However, for neutrinos the W has an axial as well as 
a vector component. Vector and axial vector meson dominance effects are expected to mainly affect the low- a; 
region. 

• Gluon recombination which can occur in a large nucleus between partons of neighboring nucleons (Nilolaev and 
Zakharov, 1975; Mueller and Qiu, 1986). This leads to an ^-dependent depletion of \ow-x partons. 

• The EMC Effect, which denotes a suppression of the structure functions from high-A targets compared to 
deuterium in the 0.2 < a; < 0.7 range. Many theoretical explanations have been developed to explain the EMC 
effect (Geesaman, et ai, 1995), which has been most clearly observed in charged lepton scattering experiments. 
These include multiquark-clusters (Krzywicki, 1976; Pirner and Vary, 1981), dynamical-rescaling (Nachtman 
and Pirner, 1984), and nuclear-binding effects (Akulinichev et ai, 1985; Dunne and Thomas, 1986). In each 
of these models, the effect is independent of the type of boson probe, and thus is expected to appear also in 
neutrino scattering using high-A targets. 

• Fermi motion of nucleons within the nucleus, which is expected to affect nuclear targets at very high x (Bodek 
and Ritchie, 1981a; Bodek and Ritchie, 1981b). 

Figure ^ shows the A-dependence of structure function data measured in deep inelastic charged-lepton scattering. 
If this A-dependence is entirely attributed to effects within the nuclear target, as opposed to propagator effects, then 
this figure should also describe the expected nuclear dependence for neutrino scattering from a high- A target. 




FIG. 8. Fit to the NMC, E665 and SLAG nuclear data. The resulting fit is used to correct the charged-lepton deuterium 
data for comparison with the GCFR result 



B. Experimental Techniques for Measuring xFs and F2 

Various experimental techniques are used to extract F2 and xF^ from the measured number of events. This section 
considers the methods used by CCFR, CDHSW and CHARM as examples. In all cases, the assumptions concerning 
the longitudinal structure function affect the final result and must be considered in any comparisons. 



1. The CCFR Measurement of F2 and xFz 

The recent CCFR E744/E770 measurements of F2 and xF^ provide an example of a technique of extracting structure 
functions in neutrino experiments with very small systematic errors. The resulting measurements demonstrate the 
high precision which recent neutrino experiments have been able to achieve. As a result of this precision, neutrino 
experiments now quantitatively test QCD at next-to-leading order. 

The structure functions F2 and xF^ are extracted from vN and vN differential cross sections given the fluxes 
(^vi $p) determined by the method described in Sec. ||(Seligman, 1997). The number of events in each x and bin 
is related to the differential cross section by: 



N = pLNa 



x — bin J Q'^—bin 



all energies 



(Pa 
dxdQ- 



■^{E)dE 



(35) 



where p is the target density, L is the target length and Na is Avogadro's number. The structure functions are 
extracted by varying the Monte Carlo input to match the number of events in each acceptance-corrected bin. 
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The CCFR results presented here are from the final analysis of the E744/E770 data which differ from the preliminary 
results of Quintas et al. (1993). Detailed information on the improvements to the analysis which lead to differences 
from the preliminary results can be found in the articles by Seligman (1997). 



om tne preuminary results can oe round m tne articles by beligman ). 
Figures ^ and |lo[ show the CCFR measurements of F2 and xF^ as a function of 



in various bins of 



Only 



statistical errors are shown, and the systematic er rors a re approximately 2% 
indicate the QCD fits which are discussed in Sec 
the QCD expectation can be probed. 



HIE 



The solid and dashed curves on the plot 
The precision of the data is such that subtle deviations from 




FIG. 9. CCFR F2 measurement. Errors are statistical only. Lines indicate the CCFR QCD fit discussed in the text. 




FIG. 10. CCFR xFs measurement. Errors are statistical only. Lines indicate the CCFR QCD fit discussed in text. 



19 



2. F2 and xFz Extracted from the Cross Section 



The structure functions can be extracted by forming the sums and differences of the differential cross section given 
in Eq. ||. This technique was used by the CDHSW and CHARM experiments. To make the measurement, the number 
of events in each x and bin is related to the differential cross section by 



1 dPa ^ N{x,Q^,E,)atot/E, 
dxdQ^ N{E^,) dxdQ^ ' 



(36) 



where N{x,Q'^ , E^) is the number of events in each (^x,Q^,E^^ bin and N{E^) is the number of events in each E^, 
bin integrated over x and . The structure functions can then be expressed in terms of the sums and differences of 
neutrino and antineutrino cross sections as 



1 dV" 

E^ dxdQ^ 
1 d^a-' 



1 



E^ dxdQ^ 
1 d^a^ 



E^ dxdQ'^ E^ dxdQ^ 



= 2g{y)xF3. 



where 



giy)^xF3 



(37) 



f{y,RL) = {l-y 

9{y) = y(i - 



2(l + i?L)'Q2' 



y^y_ 



(38) 



and Fi and xF-^ can be extracted by solving the coupled equations. 

The difference between xF^ in neutrino and antineutrino scattering, A (x-Fa), is related to the difference between 
the strange and charm sea distributions, (s — c), as shown in Eq. and discussed in Sec. 



3. Assumptions concerning Rl 



The ratio of the longitudinal-to-transverse absorption cross section, R^, can be measured from the j;-dependence 
of deep inelastic scattering data. Fits to the special function F defined as follows can be used to determine Rl. 



F{x,Q\e) 



y^GlME 



d^a'' 



dxdy dxdy 



2xFi(x,Q2) [l + eR{x,Q^)] 



(39) 



In this equation, e ~ 2(1 — y)/(l + (1 — y) ) is the polarization of virtual W boson. This equation assumes xF^ — xF^ 
and, therefore, a correction for A (xF^) must be applied. The values of Rj^ = <^lI<^t are extracted from linear fits to 
F versus e at fixed x and bins. 

There are very few measurements of i?L from deep inelastic neutrino scattering. Two experiments with published 
resuhs are CHARM (Bergsma et ai, 1984) and CDHSW (Abramowicz et ai, 1983). A comparison of the CDHSW 
measurements of Rl is made to the charged-lepton scattering results in Fig. |ri|. The best fit to the world's data on 
Rl is given by the parameterization ^^Rworid" (Whitlow et al, 1992), which is also shown in the plot. 
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FIG. 11. Measurements of i?L by neutrino and charged-lepton experiments. Solid line indicates the "worid" parameterization 
of Rl. 

CCFR E744/E770 is in the process of measuring Rl as a function of x and . The wide range of CCFR beam 
energies produces an extended range of kinematically accessible y values in each and x bin (Yang et ai, 1996). 
Therefore the lever-arm for the fit in most kinematic bins is large. 

For the F2 and xF^ results compared in the next section, various assumptions concerning Rl were made: BEBC 
and CHARM assumed Rl = 0, CCFR E616/E701 used Rl,qcd, CCFR E744/E770 implemented the Rworid param- 
eterization, and CDHSW used Rl = 0.1. 



C. Structure Functions from Neutrino Experiments 



This section compares results from precision neutrino structure function measurements. Small differences in the 
technique and assumptions used for extracting the structure functions complicate this comparison. The assumptions 
concerning the structure function Rl were discussed in the previous section. CDHSW and CCFR assume nonzero 
values for the strange sea, while CHARM and BEBC assume the strange sea to be zero. All experiments apply 
radiative corrections, but the model and order of the calculation varies in each case (De Rujula et ai, 1979; Bardin 
and Dokuchaeva, 1986). Finally, CHARM and BEBC have applied corrections at high x for Fermi smearing, while 
CDHSW and CCFR have not. The discussion below is divided between measurements using high-A tarasts and 
those from deuterium. The recent data sets presented below are available in tabular format via the Internetcl. Older 
experimental data are summarized in the tables of Diemoz et al. (1986). 



^The CCFR Collaboration E744/770 structure 
man@nevisl.columbia.edu. Data from BEBC, CCFR, 



functions are available in tabular form by contacting selig- 

, CDHSW and many charged-lepton experim ents are available in tables 

from the Durham Structure Function Archive, http: / /cpt 1 .dur.ac.uk/HEPD ATA /structurel .html| . 
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1. Heavy Target Experiments 



The experiments with high-v4 targets have the advantage of high statistics. In the comparisons between experiments 
using high-A targets, various targets were used and no nuclear corrections have been applied. 

Figure |l2| compares the structure function xF^ from experiments with iron targets as a function of Q'^ for various x 
bins. These are from CCFR E616/E701 (Oltman et ai, 1992), CCFR E744/E770 (Seligman, 1997), CDHSW (Berge 
et ai, 1991), and WHPFOR (Heagy et at, 1981). Of these experiments, CCFR E744/E770 has the smallest statistical 
error. The F2 comparison as a function of Q^, shown in Fig, p^ , indicates good agreement between the different data 
sets, except for the case of CDHSW. As shown in Sec. IIIE , the shape of the dependence as a function of x for 
the CDHSW data cannot be described by a QCD-based parameterization. The source of the disagreement is not 
understood. 
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FIG. 12. Comparison of xFs data from experiments with iron targets. 
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FIG. 13. Comparison of F2 data from experiments with iron targets. 
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FIG. 14. Comparison of xFa data from experiments with various high-A targets. 
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FIG. 15. Comparison of F2 data from experiments with various high-A targets. 

Figures |l^ and |l^ compare xF^ and F2 measured on various targets. The CCFR E744/E770 data are used to 
represent iron targets. The BEBC (Varvell et ai, 1987) data were taken on neon target. Gargamelle (Morfin et 
ai, 1981), an earlier and less precise bubble chamber experiment, used a CsHg-CFsBr mixture. The CHARM data 
were taken on marble (CaC02). The statistical errors are sufficiently large that the data are consistent, despite the 
apparent wide spread between experiments. 

The qualitative behavior of the structure functions averaged over as a function of x, shown in Fig. |l^, can be 
interpreted within QCD. Quantitative QCD studies, interpretations, and measurements will be discussed in Sec. [HE. 
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FIG. 16. Average 2:7^3 and F2 as a function of x from BEBC (open circles), CCFR E616/E701 (squares), CCFR E744/E770 
(closed circles) and WHPFOR (diamonds). 

The structure function xF^ is related to the sum of the valence quarks, and the measured distribution is peaked at 
moderate x, going to zero as x goes to zero or one. If the proton consisted of three non-interacting valence quarks, 
then this distribution would simply reflect the Fermi motion of the quarks inside the nucleon and would peak at 
X = However, because these quarks are interacting via gluons, which carry some of the fractional momentum of 
the proton, the distribution is smeared further and peaks at smaller x. The position of the peak will move toward 
smaller x with increasing due to QCD evolution. Since each neutrino experiment covers a different range of Q^, 
the data points on Fig. |l^ are not expected to agree. 

The F2 structure function is related to the sum of the quark and antiquark distributions. The F2 distribution goes 
to zero at high x and, but is approximately constant in the low- a; region. Gluon radiation and gluon splitting into 
quark- antiquark pairs results in a large number of low momentum partons that dominate the low- a; region. In fact, 
the charged-lepton scattering results on F2 from HERA have shown that at very low x, outside of the kinematic range 
accessible to present neutrino experiments, F2 rises dramatically (Aid et al, 1996b; Derrick et al, 1996a). 



2. Deuterium and Hydrogen Targets 

The BEBC experiment has used neutrino and antineutrino scattering from deuterium to obtain the best measure- 
ment of the structure functions on the proton and the neutron separately (AUasia et al., 1985; Jones et al., 1994). 
A neutrino or antineutrino interaction was identified as coming from a neutron if it had either an even number of 
prongs, or an odd number of prongs with a proton with momentum less than 150 MeV. All remaining events, with 
an odd number of prongs and hence a net total charge, were classified as interactions with protons. Misidentifications 
were corrected on a statistical basis using a Monte Carlo simulation. 
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Measurements of Neutron ard Proton by BEBC 
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FIG. 17. F^" and F^" as measured by BEBC 



The resulting structure functions for neutrino scattering on protons and neutrons, and ^2""; shown in 
Fig. 0. The most evident feature is that Fj"" ~ 2F2^ over most of the kinematic region. This is because the 
emitted in a neutrino interaction must interact with a negatively charged quark, which at high a; has the highest 
probability of being the valance d quark. Since the neutron has twice the number of valance d quarks as the proton, 
the neutron structure function is larger. These data clearly indicate the flavor sensitivity of neutrino scattering. 
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D. Comparisons of Neutrino to Charged-Lepton Experiments 
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FIG. 18. The kinematic region accessible to various charged-lepton experiments compared to the CCFR experiment 



Neutrino deep inelastic scattering experiments now have the statistical accuracy of charged-lepton scattering ex- 
periments. In this section, neutrino results are compared to charged-lepton results from SLAC (eTV), NMC {^N), 



BCDMS i^iN) and E665 (/xiV). The kinematic ranges covered by these experiments is shown in Fig. 
CCFR region shown for comparison. The HERA ep data are mainly at high and very low x, wit 
kinematic overlap with the fixed-target neutrino data. 



18|, with the 
1 only small 



1. Comparison of F2 



The parity-conserving structure function F2 is measured by both charged- and neutr al-le pton scattering experiments. 
The definitions of F2 in neutrino and charged-lepton scattering were discussed in Sec. I C. As a result of the difference 
in the couplings, a conversion must be applied in order to compare muon and neutrino experiments. To lowest order 
with the charm sea set to zero, the correction is 
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N 



1 



3(3 + s) 

5(9 + g) 



(40) 



This relation is true to all orders in the DIS renormalization scheme, which maintains Eq. |l6| as the definition of 

Assuming that the strange sea is small at high x, one expects F^^/F^^ ~^ TS ~ 0-2778 as a; — > 1. Comparing the 
CCFR data and the CDHSW data (Mishra and Sciulh, 1989) to the results from BCDMS (Benvenuti et ai, 1987) at 
high X yields a ratio of 0.278 ± 0.010 which is in excellent agreement with the expected quark charges. 

In the comparisons given below, all corrections for the differences between the charged-lepton and neutrino data 
have been applied to the charged-lepton data. The strange sea (s, s) appears explicitly in Eq. 40 and, therefore, a 
measurement of this distribution is needed for the comparison. The strange sea can be measured from the neutrino 
charm production cross section which is related to charged-current dimuon production (Bazarko et ai, 1995) described 
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in detail in Sec. IV . Add itional nuclear corrections must also be made to the charged-lepton measurements of F2, as 
discussed in Sec. " [niA2|b| . 

Figures P^P^ compare the CCFR F2 measurement to the measurements from NMC (Amaudruz et al, 1992b) and 
E665 (Adams et al, 1995) over a range of x values. The corrections for charge coupling, using the DIS scheme, 
and nuclear effects have been applied to the muon data. At low x there appears to be a disagreement between the 
charged-lepton scattering data and the neutrino data. The discrepancy decreases with increasing x until x « 0.070 
where the data are in good agreement to the highest x bins. 
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FIG. 19. Comparison of the CCFR (y) measurement of F2 to the NMC and E665 (pt) results for x = 0.0075, 0.0125 and 
0.0175. 
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FIG. 20. Comparison of the CCFR (u) measurement of F2 to the NMC and E665 (fi) results for x = 0.0250, 0.035 and 0.050 
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FIG. 21. Comparison of the CCFR (u) measurement of F2 to the NMC and E665 (fj.) results and SLAC (e) results for 
X = 0.070,0.090 and 0.110. 
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FIG. 22. Comparison of the CCFR (i^) measurement of F2 to the NMC and E665 (/.i) results for x = 0.140, 0.180 and 0.225 
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FIG. 23. Comparison of the CCFR (ly) measurement of F2 to the NMC and E665 (/i) results for x = 0.275, 0.350 and 0.450 
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FIG. 24. Comparison of the CCFR (z/) measurement of F2 to the NMC and E665 (fi) results for x = 0.550, 0.650 and 0.750. 

The high statistics data in the region of the discrepancy come from CCFR and NMC. It is possible that one or 
both experiments have underestimated their systematic errors. The most important systematic error in the CCFR 
analysis is due to the calorimeter calibration. Extensive studies of the test-beam data were undertaken in order to 
determine the muon and hadron energy calibration and the systematic error assigned to this source has been carefully 
estimated by the experiment. 

Several publications have suggested that the discrepancy is due to an incorrect strange-sea correction (Botts et 
ai, 1993; Barone et ai, 1994; Brodsky and Ma, 1996). The required distributions to eliminate the discrepancy, 
as calculated at next-to- leading order by the CTEQ collaboration (Botts et al., 1993), are inconsistent with and 
approximately a factor of two larger at low x than the direct CCFR strange-sea measurement from dimuon data 
(Bazarko et al., 1995). The CCFR measurement of the strange sea, which is to next-to-leading order and includes 
corrections for the charm-mass threshold, would have to be incorrect by 5a to account for the discrepancy. The CCFR 



strange-sea measurement and its application to the correction presented here is described in detail in Sec. IV. 

Another possible cause for the disagreement between the neutrino and charged lepton measurements of F2 may be 
the nuclear correction applied to the charged-lepton deuterium data in order to allow comparison with the neutrino 



results. As discussed in Sec. |III A 2|b| , the corrected charged-lepton data will only agree with the neutrino result if 
the entire effect is ascribed to the target material. Propagator effects, such as fluctuation of the virtual boson to 
mesons, will not be the same for neutrino and charged lepton scattering because of the additional axial component 
of the W. Such "nuclear shadowing" would be expected to affect the low- a; (low-Q^) region most strongly, thereby 
accounting for some of the disagreement. The data on these nuclear effects in neutrino and charged-lepton scattering 
are discussed in the next section. 
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2. Nuclear Effects 
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FIG. 25. Comparison of the ratio of cross sections for high-A targets to deuterium. Solid circle: BEBC data from neutrino 
- neon scattering. Open boxes: NMC data from muon - carbon scattering. 

Very little data on shadowing in neutrino scattering has been obtained. The best measurements come from the 
BEBC collaboration ratios of neon {A = 20) to deuterium (AUport et ai, 1989). This can be compared to the muon 
scattering results from carbon {A = 12) as measured by the NMC experiment (Amaudruz et at., 1995). Figure |25|a 
shows Ne/D2 directly compared to C/D2, while Fig. p5| b applies a correction factor for the theoretically expected v4~3 
dependence of the nuclear cross section to the carbon ratio (Geesaman, et al., 1995), although it should be noted that 
some experimental indicate less A dependence (Adams et ai, 1995). Data are shown for > 1 GeV^, except for the 
three lowest x BEBC data points which are from < o.2 GeV^, (0.2 < < 0.5)GcV2 and (0.5 < < l.O) GeV^, 
respectively. Within the combined statistical and systematic errors, the results from charged-lepton and neutrino 
scattering appear to agree. These results are also in agreement with an earlier, lower statistics BEBC measurement 
of this ratio which used a 3 m^ track-sensitive target filled with deuterium surrounded by neon (Cooper et ai, 1984). 
CDHSW has measured the ratio of structure functions in iron to hydrogen (Abramowicz et ai, 1984) using a 35 
m^ hydrogen tank located in front of the CDHSW detector. This experiment obtained 4,457 i^p events and 4,178 
Vp events which were compared to 50,000 neutrino and 150,000 antineutrino events in iron. Although the combined 
systematic and statistical errors are too large to establish A-dcpendcnce, the behavior of the data is consistent with 
the behavior observed in charged-lepton scattering. 

Although the NMC and E665 shadowing data which were used to obtain the nuclear correction are in statistical 
agreement, there appears to be a systematic shift between the two data sets as can be seen in Fig. ^. The nuclear 
correction is dominated by the NMC result, which has much smaller errors than the E665 data. However, if only 
the E665 data are used in a fit which is then applied to the charged-lepton deuterium results, then the CCFR/NMC 
discrepancy is reduced by approximately a factor of two for x < 0.03. However, it seems unlikely that the entire 
versus jiN low- a; discrepancy can be due to nuclear effects alone. 



E. QCD Analyses of the Structure Function Results 

An analysis of the structure function data within the QCD framework allows precise tests of the predicted QCD 
evolution and the Gross-Llewellyn Smith and the Adler sum rules. From these studies, measurements of as and the 
parton distributions can be extracted. 



1. The Strong Coupling Measured from QCD Evolution 

Perturbative QCD can predict the evolution of the structure functions given a starting set of x-dependent 
PDFs at a reference scale Qq (Altarel li and P arisi, 1977; Dokshitser, et al., 1978; Dokshitser, et ai, 1980; Gribov and 
Lipatov, 1972), as described in Sec. |HI A 1 . Therefore, an analysis of scaling violations in structure function data 



permit a determination of the parton distributions and the QCD parameter A. The CCFR E744/E770 analysis, which 
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is the most precise to date, provides an example of the extraction of these distributions and the QCD parameter at 
next-to-leadine order for four flavors in the MS renormahzation scheme (A-^^'^"'*''^). 

^ MS ' 

Two independent QCD evolution analyses were performed by CCFR. The first used only the parity violating 
structure function xF^ . The advantage of this analysis is that the evolution of xF^ is independent of the gluon 
distribution, reducing the number of free parameters. The second analysis used both xF^ and F2 data in the QCD 
fit. This increases the statistical power of the fit, but introduces extra parameters to describe the gluon distribution. 
On the other hand, the combined F2/XF3 data can constrain the gluon distribution G{x,Q'^) at moderate values of 
X. Note that charged-lepton scattering experiments cannot separate gluon distribution effects from the running of as 
as cleanly because they do not have access to the parity- violating structure function, xF^. 

The QCD fits require initial parameterizations of the parton distributions at some Qq. CCFR E744/E770 used 

xqNs{x,Ql) ^ Ansx'" - 
xqsix, Ql) = XQNsix, Ql) + As{l - (41) 
xG{x,Ql)^AGil-xy^, 

where NS, S, G refer to non-singlet, singlet, and gluon distributions, respectively; and Ajvs, As, Aq, ?7i, ?72, and 
•qo are free parameters in the fit. More complicated parameterizations of the parton distributions were found to 
yield fit results consistent with the above parameterizations (Seligman, 1997). The resulting parton distributions are 



compared to other measurements in Sec. Ill E 2 



The CCFR analysis compared data to a theoretical prediction based on the next-to-leading order QCD evolution 
program of Duke and Owens (Duke and Owens, 1984) using a technique that included the effects of target-mass, 
higher-twist, and Rl.qcd- Only data satisfying > 5 GeV^, > 10 GeV^, and x < 0.7 were included, and 
(x, Q^) bins with statistical error greater than 50% were eliminated. The effects of systematic errors were studied by 
implementing positive and negative shifts of F2 and xF^ caused by each systematic error and re-extracting A as well 
as the other fit parameters. 



The results of NLO QCD fits in the MS scheme to the CCFR structure function data are listed in Table [11 and 
are shown by the solid line on Fig. ^ and Fig. 10. The dashed lines on these figures represents an extrapolation of 



the parameterization from the fit to lower values. As expected, the fit which uses both the xF^ and F2 data has 
smaller statistical and systematic errors because of the increased number of data points used. The good per degree 
of freedom indicates consistency of the data with QCD. 

Experiment Structure Function(s) Fit A^^^*^ 

CDHSW 200+fo^ 

F2 300 ± 80 

CHARM xFs 310 ± 150 

CCFR 616/701 F2 340 ± 110 

CCFR 744/770 xFs 381 ± 55 

F2,xF3 337 ±31 

TABLE III. Values of h!l^~'^^ determined from fits to structure functions measured in various neutrino experiments. 
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The CDHSW data are not in agreement with the other neutrino measurements (see Sec. Ill CI) and with the 
expected QCD behavior. Another way to see this is to consider the logarithmic derivative of the structure functions 
with respect to as a function of x. QCD predicts the form of this distribution up to one free parameter, A, given 



the parton distributions as inputs. Figure 26 shows the CCFR logarithmic slopes and the corresponding QCD fit to 



the data along with the CDHSW data for comparison. 
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FIG. 26. CCFR and CDHSW logarithmic slopes of the structure functions. The line indicates the CCFR QCD fit result. 



An alternative and more powerful way to incorporate systematic uncertainties is to include their effects directly into 
the QCD fit. For this type of analysis, a fit to the theoretical prediction for the structure functions is compared 
to the data in each x and bin. The prediction should be compared to the data using a that includes both 
statistical and systematic uncertainties with correlations. The systematic uncertainties can be handled by introducing 
a parameter (5fc for each uncertainty k into the . Defining the structure-function vector F — (Fi^xF^) and the 
structure-function error matrix V — {(Jij), for i,j ~ F2,xF^, then the difference between the theoretical prediction 
and data and the are: 



X 



yidata 



■^theory _|_ ^ ^ (5^(F^ ■pdata'^ 



(42) 

(43) 



where F*^ is the value of the extracted structure functions when the kth systematic parameter is shifted by la. For 
the CCFR data, the results of this fitting method arc compared to the "basic" method described above in Table 
The "global fit" gives a more precise measure of A because of the constraint of QCD. 





from xF'i only 




from xFz and F2 




Fit Type: 


A± stat ± sys 


X'/DOF 


A± stat ± sys 


X'/DOF 


Basic 


387 ± 42 ± 94 


81/82 


381 ± 23 ± 58 


191/164 


Global 


381 ±41 ±37 


69/82 


337 ± 21 ± 23 


171/164 



TABLE IV. Values of hZ^-^' determined from fits to the CCFR data. The systematic error includes the experimental, 
physics model, and higher-twist uncertainties. 
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The value of A = 337 ± 31 MeV from the CCFR global fit using F2 and xF^ is equivalent to an 

a^(M|) = 0.119 ± 0.002(ea;p) ± 0.004(scak) (CCFR). 

The scale error is due to the renormalization and factorization scale uncertainties (Virchaux and Milsztajn, 1992). 
This result can be compared to other measurements as reported in Barnett et al. (1996). The CCFR result is higher 
than the previous measurement of CCFR, asiM"^) = 0.111 ± 0.005 (Quintas et al., 1993) and is also higher than 
the previous deep inelastic scattering average, as{M^) = 0.112 ± 0.005. Including the new CCFR measurement, the 
average for deep inelastic scattering experiments is 

a,(M|) = 0.117 ± 0.004 (DIS Average) 

which is lower than the LEP measurement from event shapes of 0.122 ± 0.007. However, all of the above results 
are consistent within the errors. Fig. |2^ compares the CCFR result to results from other non-neutrino experiments 
(Barnett et al, 1996). 
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FIG. 27. Measurement of as from various experiments. Results from other experiments are from Barnett et al. (1996). 



2. Parton Distributions Measured from QCD Evolution 



This section considers measurements of the parton distributions from the QCD evolution fits. The parton distri- 
butions are extracted simultaneously with the fit for A. This method can accurately determine the distributions for 
light quarks and gluons. The best method for determining the strange and charm distributions is through charged- 



and neutral-current charm production, respectively, as discussed in Sec. [V. 

The parameterizations of the parton distributions found to fit best the CCFR E744/E770 data are compared with 
several popular PDF sets available from PDFLIB (Plothow-Besch et al.^ 1993) in Fig. |2^. These PDFs are obtained 
through global fits of the world's data, including charged and neutral lepton scattering, Drell-Yan production, and 
direct photon production. All three PDF [CTEQ4M (Lai et al, 1996), MRS-R2 (Martin et al, 1996), and GRV- 
94 (Gliick et al, 1995)] have used preliminary CCFR E744/E770 structure function data in the global fits, so the 
comparison contains some correlation. The data from other neutrino experiments have not been included in most 
global fits because of poor precision or disagreement with QCD. 
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FIG. 28. The CCFR parton distributions compared to CTEQ 4M, MRS R2, and GRV 94 at = 5 GeV^ a) 
xqvaiance = xuy + xdy , b) Total Sea = ^ {xqsea + a;9sea)- 

u,d,s 

The gluon distribution is the most poorly constrained of the parton distribution measurements. This is because in 
most processes where precision measurements can be used to extract PDFs, the gluons enter only indirectly through 
the QCD predictions for the evolution. Fits to xF-^ and F2, described previously, allow one to extract the form of 
the gluon distribution. The global fit to the CCFR E744/E770 data yields the distribution 

xG{x, = 5 GeV^) = (2.22 ± 0.25)(1 - a;)4-65±o.69 (cCFR). 

Figure |2| shows this distribution as a function of x, evolved to — 32 GeV^. The shaded region indicates the 
CCFR ±1(7 errors. For comparison, the crosses show the gluon distribution as measured from jet production at the 
HI experiment (Aid et ai, 1995). The hatched region indicates the E665 gluon distribution measured from the energy 
distribution of hadrons produced in deep inelastic muon scattering (Adams et ai, 1996). Also shown are the gluon 
distributions from fits to data from a wide range of high energy experiments which are available as parton distribution 
bmctions: CTEQ4M (Lai et at, 1996) (dotted), GRV 94 HO (Gliick et at, 1995) (solid), and MRS R2 (Martin et 
ai, 1996) (dashed). The gluon distributions from the various experiments agree in the region of the CCFR data with 
X > 0.01. 
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FIG. 29. The CCFR gluon distribution evolved to = 32 GeV^ and compared to gluon distributions measured in other 
processes. Shaded region shows ilcr errors. Crosses show results from HERA jet measurements. Dotted line is CTEQ4M. 
Solid line is GRV94 HO. Dashed line is MRS R2. Hatched region is gluon distribution from E665 hadronic energy distributions. 



3. The QCD Sum Rules 

QCD predictions exist for a wide variety of "sum rules," which are integrals over quark densities. Two QCD sum 
rules can only be measured in neutrino scattering: the Gross-Llewellyn Smith (GLS) sum rule and the Adler sum 
rule. The values for these sum rules are fundamental predictions of QCD. 

a. The GLS Sum Rule. The GLS sum rule gives the QCD expectation for the integral of the valance quark 
densities. To leading order in perturbative QCD, the integral / ^{xF^) is the number of valence quarks in the proton 
and should equal three (Gross and Llewellyn Smith, 1969). QCD corrections to this integral result in a dependence 
on as (Larin and Vermarseren, 1991): 

f\Fs{x,Q')- = 
Jo ^ 

3il-^-a{n,)i^r-binf){^n (44) 

TT TT TT 

In this equation, a and b are known functions of the number of quark flavors Uf which contribute to scattering at a 
given X and Q^. This is one of the few QCD predictions that are available to order a^. 

One can determine as{Q^) from J F^dx by inverting the GLS equation. However, for a given value, there is 
a limited region in x that is accessible by the data from any one experiment. Two analysis methods address this 
problem. 

The first method extrapolates all data to a single value, resulting in data over a sufficiently wide range of x 
to evaluate accurately the integral in Eq. However, this introduces systematic errors involved in the assumptions 
for the extrapolation. To evaluate this error, three extrapolation methods can be considered: QCD evolution, l/Q^ , 
and a simple power law. The results of this analysis for several neutrino experiments are shown in Fig. The world 
average is 

j Fjdx = 2.64 ±0.06. (45) 
which is consistent with the NNLO evaluation of Eq. ^ with A = 250 ± 50 MeV. 
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FIG. 30. Summary of measurements of the GLS sum rule 



The second method evaluates the integral by combining the xF^ measurements of experiments covering different 
kinematic regions (Harris et ai, 1995). At high x, F2 is nearly equal to xF^. Therefore to improve further the 
kinematic range at high- a;, one can use F2 data from charged- lepton scattering at for x > 0.5. 

An issue which can be best addressed using the second analysis technique is that of higher-twist contributions 
to the GLS integral. Because CCFR evaluates this integral in the low-Q^ region where higher-twist contributions 
can contribute, an additional term AHT must be included in Eq. E4| (Larin and Vermarseren, 1991). This higher- 
twist contribution has been computed to be approximately 0.27/(5^3ased on three independent calculations of the 
dynamical power corrections (Braun and Kolesnichenko, 1987). The three calculations are consistent with each other 
to ~ 50%, so an error of ±0.14/Q^ should be assumed for this result. A more recent result (Ross et ai, 1994) improved 
on the technique of Braun and Kolesnichenko (1987) by explicitly retaining extra terms dropped in the earlier analysis. 
As a result, the extracted value is more reliable, leading to a smaller theoretical error, (0.16 ± Q.01)/Q^ (Ross et ai, 



b. The Adler Sum Rule. The Adler Sum Rule predicts the difference between the quark densities of the neutron 
and the proton, integrated over x. To leading order: 



The only assumption introduced in this calculation is that the quark and antiquark seas are equal; there is no 
assumption concerning the equality of, for example, the u and d seas. 

Using the F2^ and F^" data shown in Fig. |l^, BEBC has demonstrated that the Adler sum rule is consistent with 
2.0 (AUasia et ai, 1985). For this review, we have evolved the Fj"" — measurements to = 5 GeV^ as shown 
in Fig. 1^ (closed points). We then evaluate the Adler sum rule under two assumptions for the low-x behavior. For 
the first method, we assume that the lowest bin is a constant equal to the value of the lowest-x data point. In this 
case, the integral is 2.05 ± 0.15 at = 5 . Alternatively, for the second method, we fit for the low- a; behavior using 
the polynomial shown by the solid line. Using this polynomial to evaluate the integral for < a; < 0.01, we obtain 
1.87 ± 0.15, which is also within Icr of the expectation. 



1994). 




(46) 
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FIG. 31. Evaluation of the Adler Sum Rule using the data from BEBC evolved to = 5 GeV^ (closed circles). The solid 
line shows the fit to the low-a; data. The open circles show the integral from each bin to a; = 1. Using the fit to evaluate the 
integral for x < 0.01, the sum rule is found to be 1.87 ± 15. 

The comparable sum rule for charged-lepton scattering is the Gottfried Sum Rule (Gottfried, 1967): 

[\Fr{x,Q')-F^^{x,Q'))- = 
Jo X 

^ {u{x) - d{x)){el - el)dx = i (47) 

where and are the charges of the u and d quarks. Inherent in this form of the Gottfried sum ndc is the 
assumption that the u distribution is equal to the d distribution. There is no a priori reason to expect a difference in 
these distributions and there was therefore some surprise when the NMC experiment measured the Gottfried sum to 
be 0.240 ± 0.016 (AUasia et at, 1990). This result could be explained ifu^d. Note that this assumption u = d was 
not required in the Adler sum rule. 



IV. NEUTRINO CHARM PRODUCTION: MEASUREMENTS OF THE STRANGE SEA, CHARM SEA, 

AND \Vcd\ 

A. Introduction 



As shown in the previous sections, neutrino DIS is particularly well suited for measuring the parton densities 
due to the neutrino's ability to resolve the flavor of the nucleon constituents. In addition, neutrino scattering is an 
effective way to study the dynamics of heavy quark production, due to the light to heavy quark transition at the 
charged current vertex. In particular, neutrino charm production can be used to isolate the nucleon strange quark 
distributions, xs{x) and xs{x), and study the transition to the heavy charm quark. The strange quark distribution 
function is of particular theoretical interest since it may contribute to the low-Q^ properties of the nucleon in the 
non-perturbative regime. Charm production is also an important testing ground for NLO perturbative QCD due to 
the large contribution from gluon initiated diagrams. In addition, understanding the threshold behavior associated 
with the heavy charm mass is critical to the extraction of the weak mixing angle, sin^ 6w, from neutrino neutral 
current data as described in the next section. 

One distinctive signature for the production of charmed quarks in neutrino- and antineutrino-nucleon scattering 
is the presence of two oppositely-charged muons. (Studies of same-sign dimuon production have shown that the 
observed signal is completely dominated by non-prompt background sources which preclude any quantitative physics 
measurements (Sandler et ai, 1993).) In the case of neutrino scattering, the underlying process is a neutrino interacting 
with an s or d quark, producing a charm quark that fragments into a charmed hadron. The charmed hadron's 
semileptonic decay {BR =?« 10% ) produces a second muon of opposite sign from the first. 

i/^ N — ^ n'+c + X (48) 
s -h Ai+ -I- i/^ 
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The analogous process with an incident antineutrino proceeds through an interaction with an s or d antiquark, 
again leading to oppositely-signed muons in the final state. 

+ N ^ ^++c + X (49) 

The heavy charm quark is expected to introduce an energy threshold in the dimuon production rate. This effect 
has been described in the past through the "slow rescaling" model (Barnett, 1976; Georgi and Politzer, 1976), in 
which ^"^j the momentum fraction carried by the struck quark, is related to the kinematic variable x = Q'^ /2Mv 
by the expression = x{l + m?^/Q^) where is the mass of the charm quark. (A more complete treatment 
of the kinematics associated with heavy quark production modifies this model slightly, as discussed below.) Charm 
production by neutrinos (or anti-neutrinos) is a substantial fraction of the total cross section and grows from threshold 
to about 10% at high energy as shown in Figure p^. 
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FIG. 32. Fraction of the neutrino (anti-neutrino) total cross section that is associated with charm production. The dashed 
curves are from a leading-order calculation with rric = 1.31 GeV and the solid are from the same calculation with = 0. 



Neutral current scattering off the charm sea can be isolated and studied by looking for the production of wrong-sign 
single muons (WSM) in the interaction of a relatively pure beam. This is one of the only techniques for probing the 
charm content in the nucleon which is assumed in the standard calculation to evolve from zero at a threshold ~ 4 
GeV^ . The process for study is neutrino neutral current scattering off a charm quark followed by the semi-leptonic 
decay of the charm quark into a wrong-sign muon. The size of this process is therefore dependent on the number and 
distribution of charm quarks in the nucleon. The main backgrounds come from the contamination in the beam 
and decays of secondary pions or kaons in the hadronic shower. The signal is smaller than the background and at 
present only upper limits on the size of the charm sea are available. 



B. Differential Cross Section for Dimuon Production and Heavy Quark Effects 

The differential cross section for dimuon production is expressed generally as 

dM^.^^.-.-X) ^ dM^N-. cX) ^^^^ ^^^^ ^ ^^^^ 
at, dy dz dt, dy 

where the function D(z) describes the hadronization of charmed quarks and is the weighted average of the semi- 
leptonic branching ratios of the charmed hadrons produced in neutrino interactions. 

The heavy charm quark introduces an energy dependence in the charm production rate and also affects the measured 
angular (or y) distribution compared to that expected for a massless quark. (Additional thresholds related to the 
masses of final-state charmed mesons and baryons are present but only become important near production threshold.) 
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These are kinematic effects for which ^, the momentum fraction of the struck quark, is related to the Bjorken scaHng 
variable x, through the expression (Aivazis et ai, 1994) 




m^ + mg + A(- Q^m^,mg) 

, (51) 



2g2 

where rric is the charm quark mass and refers to the initial state quark mass, either the strange quark or the down 
quark, and A(— Q^, ml, ml) is the triangle function, defined by A(a, 6, c) = a?- + IP- + — 2{ab + bc + ca). The full 
expression for ^ can be simplified by neglecting the small effect of the initial state quark mass to yield 

e.J.,'4)U-^^)^ (52, 



Relating ^ and x through the charm quark mass is referred to as "slow-rescaling" (Georgi and Politzer, 1976; Barnett, 
1976). A comparison of the full ^ calculation to the simple ^° expression given above yields a difference of 10% (2%) 
for X = 0.2 and — 1.0 (5.0) GeV^; this difference is much smaller at lower x and higher Q^. 

At leading order charm is produced by scattering directly off of strange and down quarks in the nucleon. The LO 
differential cross section for an isoscalar target, neglecting target mass effects, is given by: 

(Pa(iy^,N cX)] GlME^ , _ 2n . 2m ,t. ,2 
{ =^TTr75277W{Kw^,M' 

d^dy ) LO T^i^ + Q /M^r 



+ 2^.(^,M^)l^.in^l-^j, (53) 

where ^u{^,fjp), ^c?(^,/i^) and ^s(C,^^) represent the momentum distributions of the u, d and s quarks within the 
proton (the corresponding 77^ process has the quarks replaced by their antiquark partners). \Vcd\ and \Vcs\ are the 
CKM matrix elements. The dependence of the parton distributions on the scale /i^ is specified by QCD (Altarelli and 
Parisi, 1977; Dokshitser, et al., 1978; Dokshitser, et al., 1980; Gribov and Lipatov, 1972), but the precise dependence 
of iJp on and x is somewhat arbitrary in finite-order perturbation theory. The finite charm mass included in the 
(1 - ml/2MEyCl expression induces an effective Callan-Gross violation such that jiLO-charm^^^ q2^ ~ 1 + ml/Q'^. 
In the modified leading order analysis of Rabinowitz et al. (1993), Callan-Gross violation is included by replacing the 
term (l - ml/{2ME^S,)) in Eq.(|3|) with 

1 - ml/{2ME,0 - ll^2miQ? " ^ " ^^^/(^^'-^ + 

using external measurements of the structure function Rl{£,,Q^) (Dasu et ai, 1988). In the next-to-leading or- 
der (NLO) formalism, violation of the Callan-Gross relation emerges as a consequence of the more complete QCD 
calculation. 

The LO expression illustrates the sensitivity of the process to the strange quark sea. Charm (anticharm) production 
from scattering off d (d) quarks is Cabibbo suppressed. In the case of charm produced by neutrinos, approximately 
50% is due to scattering from s quarks, even though the d quark content of the proton is approximately ten times 
larger. In the case of antineutrino scattering, where d quarks from the sea contribute, roughly 90% is due to scattering 
off s quarks. 



C. Next-to-Leading Order Corrections 



Because neutrino charm production has a large sea quark component at leading-order, the next-to-leading-order 
gluon- initiated contributions are significant (Aivazis et ai, 1990). Gluon-initiated production of charm proceeds 
through both the t and u channels as shown in Figure ^3|a. The size of the gluon distribution, which is an order 
of magnitude larger than the sea quark distribution, compensates for the extra power of as involved in the gluon- 
initiated diagram. The NLO quark-initiated diagrams, shown in Fig. |33|b, in which a gluon is radiated, also enter 
the perturbative expansion at 0{as), but the corrections from these diagrams to the cross section are small. Several 
calculations including the next-to-leading-order formalism have been done (Aivazis et aL, 1994; Kramer and Lampe, 
1992; Van der Bij and van Oldenborgh, 1991). The treatment of the gluon-initiated diagrams in these has been 
cross-checked and found to be consistent. 
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FIG. 33. Mechanisms that contribute to neutrino production of charm up to 0{as)- a) The dominant diagrams: the 
leading-order quark-initiated diagram, and the t channel and u channel gluon-initiated diagrams, respectively, b) The radia- 
tive-gluon and self-energy diagrams. 

As with all applications of perturbative QCD, a theoretical uncertainty is associated with the choice of factorization 
and renormalization scales. Some scale dependence is unavoidable for any calculation done to finite order in as. Most 
analyses assume that the factorization and renormalization scales are both equal to a single parameter, /Lt. The fi 
scale is interpreted as setting the boundary between the coUinear and non-coUinear regions of the p± integration over 
the final states where p± is the transverse momentum of the initial state quark coming from the gluon splitting. 
Therefore, a scale proportional to p'^^^ = A{W^ , iir^, M^) / V4W^ is suggested by Aivazis et al. (1994). Figure pi] 
shows the scale dependence of the differential cross section for charm production, where the abscissa is in units of 
pmax_ rp-j^g scale dependence is weak for values above one unit of p™'^^, but there is a strong dependence when fi is 
below this value. 
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FIG. 34. The scale dependence of the differential cross section for neutrino and antineutrino production of charm, where 
identifies the factorization and renormalization scales. The scale on the abscissa is in units of p™^". For E = 200 GeV and 
y = 0.5, the x = 0.05, 0.15, 0.25 lines correspond to = 6.6, 6.2, 5.8 GeV/c, respectively. 



D. Experimental Issues for Dimuon Measurements 



1. Data Samples and Event Selection 



Many experiments have investigated neutrino and antineutrino production of charm since it was first observed in 
1974 (Benvenuti et at, 1975). The experiments considered here include the highest statist ics experiments using each of 



the main experimental techniques. Most of the experiments have been described in Sec. ID and only details relevant 
to charm measurements will be given here. The largest data samples are available from the CCFR (Bazarko et ai, 
1995) and CDHS (Abramowicz et ai, 1982) experiments. The CCFR collaboration used a high-energy, quadrupole- 
triplet neutrino beam and a massive, high-density steel detector at the Fermilab Tevatron. The CDHS collaboration 
at CERN used a magnetized steel target-calorimeter illuminated by a horn-focused neutrino beam with a somewhat 
lower energy spectrum. Results are also available from the Fermilab 15-ft bubble chamber experiment (15-ft BC) 
for fj,~e~^ events (Baker et ai, 1991). In this experiment, the low threshold, Pe+ > 0.3 GeV, combined with good 
statistics for < 30 GeV give good sensitivity to the slow-rescaling threshold behavior. The final data set considered 
is the E531 neutrino emulsion experiment (Ushida et ai, 1988) at Fermilab which measures inclusive charm particle 
production directly by identifying the charm decay in the emulsion. Thus, this experiment can provide information 
on the charm species produced by neutrinos and the kinematics of the production. Table ^ summarizes the neutrino 
charm production data samples from these experiments. 



Experiment 


E„ (GoV) 


/i events 


/i^ I events 


Background (%) 


CCFR (Bazarko et ai, 1995) 


30 - 600 


5030 


1060 


15 


> 5 GeV 


(> 100) 


(3721) 


(493) 




CDHS (Abramowicz et ai, 1982) 


30 - 250 


11041 


3684 


13 


> 5 GeV 


(> 100) 


(3589) 


(452) 




15-ft BC (Baker et ai, 1991) 


1 - 200 


461 




18 


p,+ > 0.3 GeV 










E531 (Ushida et ai, 1988) 


1 - 250 


122 




3 


v — Emulsion 




(Charm Events) 







TABLE V. Summary of data samples from some neutrino charm production experiments. The CDHS and CCFR experiments 
detect dimuon events, the 15-ft BC measures jj, — e events, and the E531 experiment detects the inclusive decays of charmed 
particles in an emulsion target. 
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2. Non-prompt and Other Background Sources of Dimuons 



For the CCFR and CDHS experiments, non-prompt pion and kaon decays constitute the main background to 
the charm-initiated signal. The high density of their target-calorimeters minimizes this contamination due to the 
short interaction length of the detector. The background is estimated from Monte Carlo calculations and test beam 
measurements of muon production in hadron showers. The systematic uncertainty in this background has been 
estimated (Sandler et ai, 1993) by the CCFR collaboration to be 15%(20%) for i/(l7) induced showers and mainly 
comes from uncertainties in the secondary particle production spectrum including heavy nuclear target effects. For 
the 15-ft BC experiment, the background is from Dalitz decays of vr" mesons and other 7 conversions. Finally, the 
E531 emulsion experiment has a small background from non-charm events that pass their charm identification cuts. 
The background level for each experiment is listed in Table 0. 



3. Charm quark production, fragmentation, and experimental acceptance 

The production and fragmentation of the charm quark into charmed hadrons must be modeled in order to extract 
information from the observed dilepton events. For the LO formalism, the kinematic suppression due to the mass of 
the charm quark, rric, is included in Eq. ^ and the charm quark direction is defined by the W'^ q ~^ c kinematics. 

The gluon-initiated diagrams in the NLO calculation proceed by W~^g cs as shown in Fig. |33|a. This NLO 
production proceeds through both the t and u channels which dominate different regions of phase space. In the t 
channel, the gluon splits into an ss pair and the c quark emerges from the VF-boson vertex. In the u channel, the 
legs of the c and s quarks are crossed — the gluon splits into a cc pair and the s quark emerges from the VF-boson 
vertex. In the VF-boson-gluon center of mass frame, the c quark is produced at an angle 9* relative to the VF-boson 
direction. The production angle is related to the momentum of the c quark in the lab. When 9* is small — t channel 
dominance — the c quark carries most of the W^-boson momentum. As 9* approaches tt — u channel dominance — the 
c quark emerges with little momentum in the lab. Acceptance corrections which are dependent on 9* need to be 
included in the analysis of any experimental measurement; these corrections are largest at high 9* where the decay 
products of the outgoing c quark are hard to detect. 

The charm quark fragmentation uncertainties are a second important source of systematic error for the dilepton 
physics measurements. The fragmentation of the charm quark into D mesons is typically parameterized using the 
Collins and Spiller function (Collins and Spiller, 1985), 

D{z) ^N[{1- z)/z + ec{2 - z)/{l - z)] (1 + z)'[l - (1/z) - ec/{l - zT'; (55) 

older measurements also employed the Peterson function (Peterson et ai, 1983), 

D{z) = Nz[l - (1/2) - ep/(l - z)]2"\ (56) 

where z = pd/p'd^^ the fraction of its maximum possible momentum that the D meson carries and the e parameters 
determine the shape of the distribution. In principle, e"^e~ collider data could provide precise data on the e parameters. 
However, because of uncertainties due to the different quark environment in e"'"e~ relative to neutrino scattering, the 
e parameters in the fragmentation models are best determined from fits to neutrino charm production data. Fits 
to the charm particle spectrum in the E531 data (Ushida et ai, 1988) give ep = 0.18 ± 0.06 and fits to the CCFR, 
z^is = EiJ{Ei^ + Ehad) distribution yield ec = 0.81 ± 0.14 (Bazarko et ai, 1995) (or ep = 0.22 ± 0.05 (Rabinowitz 
et al, 1993)). These can be compared to the CLEO measurements (Bortoletto et al, 1988): The Collins and Spiller 
function fits the CLEO data well with ec = 0.64±0.14, but the Peterson function does not reproduce the data for any 
value of ep . The CDHS collaboration in their dimuon analysis (Abramowicz et ai, 1982) parameterized the charm 
fragmentation phcnomcnologically by functions ranging from a delta function ai z = 0.68 to a flat distribution. 

The semileptonic decay of the charmed particles also needs to be modeled and the effect of missing energy from 
the outgoing neutrino needs to be included. For this modeling, the three-body decay kinematics is taken from e^e^ 
collider studies including the differences for the various charmed particle species. The composition of charmed particles 
produced in neutrino interactions has been measured by the E531 j/— emulsion experiment (Ushida et al, 1988). With 
an El/ > 30 GeV cut, the production is dominated by neutral and charged D mesons. From a recent re-analysis 
(Bolton, 1994a; Bolton, 1994b) of the neutrino produced charm particle data, the production fractions are 0.58 ±0.06, 
0.26 ± 0.06, 0.07 ± 0.05, and 0.07 ± 0.04 for 0°, D+, D+, and A+ respectively for this energy range. 
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E. Dimuon Measurements and Results 



1. Charm Production Rate 



The energy dependence of the neutrino charm production rates provides a direct test of the modehng of the heavy 
charm quark production and the thresholds appropriate for the charm particles in the final states. At low energy, 
El, < 20 GeV , quasi-elastic scattering, such as ti ^ /^^ A+ , is a large fraction of the charm production and 
accounts for more than 25% of the charm cross section. Above 50 GeV, the exclusive final state effects become small 
(below 5%) and the simple parton model approach in which the production and hadronization factorizes becomes 
valid. Comparing the predicted energy dependence expected from the "slow-rescaling" model with measured rates 
corrected for acceptance, smearing, and other kinematic cuts shows good agreement between the data and model for 
the CCFR, 15-ft BC, and E531 experiments as shown in Fig. The CDHS data, however, is somewhat below the 
prediction and the other data, especially at low energy. 
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FIG. 35. The dilepton to single muon rate, corrected for acceptance, smearing, and minimum second lepton cut versus 
neutrino energy for the CCFR, CDHS, and 15-ft BC experiments. Also, shown is the E531 emulsion experiment where the 
charm rate has been multiplied by a semileptonic branching ratio of 0.11. The dashed curve is a LO "slow-rescaling" model 
prediction with m^ = 1.31 GeV. (The errors on the CDHS points include uncertainties in fragmentation and normalization; the 
other points do not.) 



After correcting for the "slow-rescaling" threshold corresponding to an appropriate nic, the rates become less 
dependent on E,y, exhibiting only the sharp, low energy threshold behavior associated with the production of heavy 
charmed mesons [VF^ > {Md + Mproton)'^]- The CCFR data shown in Fig. ^ is in good agreement with this model 
for a charm quark mass of rric — 1.31 GeV/c^. In the next section, the energy dependence of charm production will 
be used to determine nic. 
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FIG. 36. Opposite-sign dimuon rates versus for the CCFR v^j, (top) and Z/j^ (bottom) data. Rates corrected for acceptance, 
smearing, and kinematic cuts are indicated by the squares. Those corrected for slow rescaUng with mc = 1.31 GeV are given 
by the circles. The curves indicate the LO model prediction with mc = 1.31 GeV before (dotted) and after (dashed) correcting 
for the finite charm mass but including charm particle mass effects. 



2. Extraction of the Strange Sea, xs{x,Q^), and Charm mass, mc- 

The CCFR collaboration have done both a leading order (Rabinowitz et ai, 1993) and next-to- leading (Bazarko et 
ai, 1995) order analysis of their diinuon data samples. For these analyses, the single-niuon and dimuon events were 
simulated using Monte Carlo techniques and compared to the observed dimuon distribution to extract the various 
physics parameters. Quark and antiquark momentum densities were obtained from the measured CCFR structure 
functions (Quintas et ai, 1993; Leung et ai, 1993) using the following procedure. 

The F2 and xF^ structure functions are used to determine the singlet, xqsi{x, jj}^) = xq{x,fi^) + xq{x,iJ.'^), non- 
singlet, xqNs{x, n^) = xq{x, ii'^) — xq{x, fj,"^), and gluon, xg{x, ji^), distributions. These distributions are obtained from 
leading-order or next-to-leading-order QCD fits to the structure function data, described in Seligman (1997), using 
the QCD evolution programs of Buras-Gaemers (Buras and Gaemers, 1978) or Duke and Owens (Duke and Owens, 
1984) respectively. 

To resolve the strange component of the quark sea, the singlet and non-singlet quark distributions are separated 

by flavor. Insofar as isospin is a good symmetry, the experiment is insensitive to the exact form of the up and down 
valence and sea quark distributions, because the neutrino target is composed of iron which is nearly isoscalar. An 
isoscalar correction accounts for the 5.67% neutron excess in the iron target. The proton valence quark content, 
xqv{x,^^) = xqNs{x, 1^"^), is parameterized by 

xqv{x, fi^) = xuv{x, fi^) + xdv{x, fi^), 

xdv{x, n^) = A4{1 — x)xuv{x, n^), (57) 

where the shape difference for xdv{x) better fits charged- lepton scattering measurements of F2/F2 (Amadruz et at, 
1992a). Ad is fixed by demanding that the ratio of the number of d to u valence quarks in the proton is 1/2. The 
non-strange quark and antiquark components of the sea are assumed to be symmetric, so that xu{x, fi^) = xus{x, fi^), 
xd{x,fi'^) = xds{x,iJ,^). The isoscalar correction is applied assuming xu{x,^'^) = xd{x,^'^). 
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The strange component of the quark sea is allowed to have a different magnitude and shape from the non-strange 
component. The strange quark content is set by the parameter 



J^[xs{x, fi'^) + xs{x,fi'^)] dx 
[xu{x, /i^) + xd{x, /i^)] dx 



(58) 



where k — 1 would indicate a flavor SU(3) symmetric sea. The shape of the strange quark distribution is related to 
that of the non-strange sea by a shape parameter a, where a = would indicate that the strange sea has the same 
X dependence as the non-strange component of the quark sea. Assuming that xs{x,^^) and xs{x,^^) are the same, 
the sea quark distributions are then parameterized by: 



xq{x, /i^) = 2 



xu{x, y?) -f xd{x, /i^) 



xs{x, fi^) = As{l — xY 



xs(x, fl^ 



xu{x, -I- xd{x, fi^) 



(59) 



where As is defined in terms of k and a. 

A y"^ minimization is performed to find the strange sea parameters k and a, the values of Be and rric that appear in 
Eq.'s 50 and [5^, and the fragmentation parameter e, by comparing the data and Monte Carlo x^is, E^is and Zvis event 
distributions. Projections of data for each of these variables are compared with the NLO fits (Bazarko et ai, 1995) in 
Figs. H and M Taking \Vcd\ = 0.221 ± 0.003 and \Vcs\ = 0.9743 ± 0.0008 (Barnett et at, 1996) as input values, 



the extracted CCFR LO and NLO parameters with their statistical and systematic errors are presented in Table Vl 





Fragmentation 


x7dof 




a 


Be 


rUc (GeV) 


CCFR (Bazarko et al, 1995) 


Collins- Spiller 


52.2/ 


0.477 


-0.02 


0.1091 


1.70 


NLO fit 


e = 0.81 ±0.14 


65 


±0.045 ± 0.024 


±0.57 ±0.27 


±0.0078 ± 0.0057 


±0.17 ±0.09 


CCFR (Bazarko et al, 1995) 


Peterson 


41.2/ 


0.468 


-0.05 


0.1047 


1.69 


NLO fit 


ep = 0.20 ±0.04 


46 


±0.053 ± 0.025 


±0.47 ±0.27 


±0.0076 ± 0.0057 


±0.16 ± 0.11 


CCFR (Rabinowitz et al, 1993) 


Peterson 


42.5/ 


0.373 


2.50 


0.1050 


1.31 


LO fit 


ep = 0.20 ± 0.04 


46 


±0.045 ± 0.018 


±0.58 ±0.31 


±0.007 ± 0.005 


±0.21 ±0.12 


CDHS (Abramowicz et al, 1982) 


(5-function 




0.478 


Not fit 


0.0839 


Not fit 


LO fit 


(z = 0.68) 




±0.094 


Assumed 0.0 


±0.0147 


Varied 1.2 - 1.8 



TABLE VL Next-to-leading-order and leading-order fit results, assuming xs{x) — xs{x). Errors are statistical and system- 
atic, except that the errors on the fragmentation parameters are statistical only. 
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37. 


The 




distribution for the CCFR v- and !7-induced dimuon events. Data are given by the points and the soUd 
histogram is the result of fitting the dimuon event simulation. The dotted histogram is the background contribution to the 
former from pion and kaon decay. 





E™ (GeV) 



FIG. 38. The Evis distribution for the CCFR v- and T/-induced dimuon events after correction for ■k/K decay backgrounds 

and V — V misidentification. Data arc given by the points and the solid histogram is the result of fitting the dimuon event 
simulation. The quarks contributing to the histogram are indicated by the curves where s, d„, and d^ are the strange sea, d 
valence, d sea components, respectively. 
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FIG. 39. The x^is distribution for the CCFR v- and i7-induced dimuon events after correction for tt/K decay backgrounds 
and V — V misidentification. Data are given by the points and the solid histogram is the result of fitting the dimuon event 
simulation. The quarks contributing to the histogram are indicated by the curves where s, d^, and ds are the strange sea, d 
valence, d sea components, respectively. 



The CDHS collaboration (Abramowicz et ai, 1982) extracts the combination \Vcd\'^Bc from the slow-rescaling and 
acceptance corrected dimuon rates at high energy, E^, > 80 GeV. The rate is used to remove the strange quark 
contribution to the i/^ rate, leaving the down to charm quark contribution as follows: 

3 

where = ^I'^v) / '^I'^u) ^nd r = gI^IgY = 0.48 ± 0.02. 

In the CDHS analysis, the combination k| Vcsp/| Vcdpis found by a shape analysis of the observed x^s distribution 
for neutrino induced dimuons where the distribution is separated into scattering off d ox s quarks. The antineutrino 
dimuon data gives the strange quark distribution and their measured charged-current structure functions gives the 
distribution for down quarks. After including the fragmentation uncertainties, a 10% scale error, and variations in 
the charm mass between 1.2 and 1.8 GeV/c^, the results are: 

I VcdpSc = (4.10 ± 0.70 ± 0.16) X 10^3 (CDHS), 
W P 

= 9.3± 1.6±0.9 (CDHS). 

where the first error is the total error except the rric variation which is given as the second error. Using the same 
CKM matrix elements as above, \Vcd\ = 0.221 ±0.003 and \Vcs\ = 0.9743 ± 0.0008 (Barnett et at, 1996), these values 



can be converted to the measurements of Br and k shown in Table VI 



3. Discussion of Results 



Parton distributions are defined to a given order and scheme in QCD. Therefore, the magnitude of a given parton 
distribution differs between leading-order and next-to-leading-order. The parameters a and rric show shifts between 
LO and NLO, but k and Be are similar. For the CCFR NLO analysis, the nucleon strange quark content is found to 
be K = 0.477 ± 0.051, indicating that the sea is not SU(3) symmetric. This is qualitatively the same result as from 
the LO analysis of CCFR and CDHS. 
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Since a nonzero value of a would indicate a shape difference between xq{x) and xs{x)^ the CCFR NLO value 
a = —0.02 1 Q indicates no shape difference at NLO. At leading order, CCFR finds the strange quarks softer than 
the overall quark sea. This comes about in part because the CCFR LO analysis includes a longitudinal component, 
Rl{£,, Q^) , in both the inclusive charged-current and charm production differential cross section as shown in Eq. |^. 
This procedure is an approximation and does not correctly treat the differences in the longitudinal component for 
heavy and light quark production. On the other hand, the NLO formalism does include these mass effects correctly 
and should give reliable measurements of the xq{x) and xs{x) seas. Fig. ^ compares the results for the quark and 
strange seas determined by the NLO CCFR analysis which shows that the magnitude is different but the shape is the 
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0.0 
0.01 
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FIG. 40. The quark and strange sea distributions, xq{x,fj,'^) and xs{x,^'^), at /x^ = 4GeV^ determined from the CCFR 
next-to- leading order analysis. The band around the NLO strange sea distribution indicates the ±la uncertainty in the 
distribution. 



The charm quark mass parame ter from the CCFR NLO fit is 1.70 ± 0.19 GeV which differs from the leading- 
order result. As shown in Section IVEl, the LO analysis with the "slow-rescaling" correction for the heavy charm 
quark is able to parameterize the threshold behavior of the measured process. On the other hand, the NLO value of 
TOc, while also describing the threshold behavior well, correctly includes the kinematic effects associated with heavy 
quark production and provides a better comparison with measurements derived from other processes involving similar 
higher-order perturbative QCD calculations. For example, the photon-gluon-fusion analysis of photoproduction data 
of Anjos et al. (1990) finds TOc — l-74^g :[g GeV which is in agreement with the NLO neutrino value. 

Gliick, Kretzer, and Reya (Gliick et al., 1996) have erroneously claimed that the acceptance correction was applied 
incorrectly in the CCFR NLO analysis (Bazarko et al., 1995). In the analysis procedure of Aivazis et al. (1994), the 
cross section is divided into three terms: LO, NLO, and subtraction. The CCFR analysis consistently corrects for 
all experimental acceptance effects using the proper kinematics for each term. For the LO and subtraction term, the 
underlying distribution corresponds to the appropriate W'^s c kinematics. The acceptance for the NLO term uses 
the W'^g cs kinematics with the calculated matrix element from Aivazis et al. (1994). Gliick, Kretzer, and Reya 
also claim that the NLO to LO difference is not supported by the currently available parton distribution functions 
(Lai et al., 1996; Gliick et al., 1995; Martin et al., 1996). As stated above, the CCFR LO analysis includes higher 
order longitudinal cross section components that make their parton distributions different from strict LO analyses. 
The most recent and complete CCFR NLO analysis (Seligman, 1997) gives parton distributions that are consistent 
with those from the various global fits as shown in Fig. |2^. It is interesting to note that the most recent GRV 94 
(Gliick et al., 1995) strange sea distributions are much smaller than the CCFR NLO measurements and the other 
global fit results as shown in Fig. Publication in the future of acceptance and smearing corrected neutrino charm 
production cross sections or structure functions by the CCFR group will allow independent analyses and global fits 
using the different NLO calculations. 
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FIG. 41. The CCFR NLO strange sea distributions compared to CTEQ 4M, MRS R2, and GRV 94 at fj.^ = 5 GeV^ 

The CCFR collaboration has also investigated the dependence of their NLO results on factorization and renormal- 
ization scale uncertainties. Their standard results quoted in Table ^ use /i — 2p™^^ for both scales and find the scale 
uncertainty by varying /i between p™^'^ and Sp^J^. Fit results with various choices of the common factorization and 



renormalization scale are presented in Table VII. These results indicate that the data favor /^^ scales with smaller 



magnitudes but that the values of the fit parameters are fairly insensitive to the choice of scale. 



choice of scale, fi^ 






a 




m, (GeV) 




50.4 


0.513 


0.18 


0.0987 


1.71 




52.2 


0.477 


-0.02 


0.1091 


1.70 




54.4 


0.460 


-0.10 


0.1142 


1.68 


Q' 


51.7 


0.423 


-0.37 


0.1074 


1.80 


(2Q)' 


56.1 


0.410 


-0.46 


0.1159 


1.71 


mf 


59.4 


0.408 


-0.54 


0.1206 


1.73 




52.5 


0.421 


-0.03 


0.1066 


1.65 




57.1 


0.409 


-0.16 


0.1154 


1.64 


+ {2mcf 


52.8 


0.428 


0.00 


0.1068 


1.62 




57.3 


0.415 


-0.15 


0.1161 


1.63 



TABLE VII. Central values of the CCFR fit parameters for various choices of the QCD scale p? . Each fit contains 65 degrees 
of freedom. 
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4- Comparisons to other strange sea measurements 



The comparison of neutrino and charged- lepton measurements of the F2 structure function (See Section HID) can 
also be used to study the magnitude of the strange sea. In the DIS renormaUzation scheme, the strange sea is related 
to the two F2 measurements by 

i {X S{X, + 2: Q2)) ^ ^F^^N q2) _ 3 p^^N q2) (g^) 

The CTEQ collaboration in their CTEQ-IMS parton distributions (Botts et ai, 1993) tried increasing the strange 
sea at low x to be consistent with Eq. The strange sea from this procedure is much larger than the direct 

CCFR NLO measurements (Bazarko et a/., 1995) and other parton distributions (Botts et ai, 1993; Martin et ai, 
1993; Gliick et ai, 1992) as shown in Fig. |4^. It therefore seems unlikely that a larger strange sea is the explanation 
for the neutrino/charged-lepton difference at low x. The same conclusion has been put forward by Gliick et al. (1996) 
who state that if the CCFR (Seligman, 1997) and NMC (Amaudruz et ai, 1992b) data are correct, this discrepancy 
"will constitute a major problem which cannot be solved within our present understanding of the so far successful 
perturbative QCD." 
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FIG. 42. Strange quark distributions xs{x,fi^) from the CCFR experiment and previous global fits by CTEQ, MRS and 
GRV for = 4 and 20 GeV^ 



5. Tests of xs(x) 7^ xs(x) 



In principle, the momentum distributions of s and s quarks need not be the same as long as the strangeness content 
of the nucleon is constrained to have equal numbers of s and s quarks. Theoretical work has explored the possibility 
that the nucleon contains a sizable heavy quark component at moderate x — the possibility of so-called intrinsic heavy 
quark states within the nucleon (Brodsky et al., 1981; Brodsky and Ma, 1996). Postulating intrinsic strange particle 
states, such as a K component of the proton wave function, leads to the prediction that the s quark momentum 
distribution will be harder than the s quark distribution (Burkardt and Warr, 1992). 

The CCFR collaboration has explored this possibility by performing an analysis (Bazarko et al., 1995) in which the 
momentum distributions of the s and s quarks are allowed to be different. For this study the sea quark distributions 
are parameterized by: 
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2\ I xu{x, + xd{x, ^^)\ xs{x, fi^) + xs{x, fj,'^) 



xq{x, fi ) — 2 
xs{x,n^) = As(l-a;)" 
xs{x,fj.^) ^A'.il-x)"' 



xu{x, fj,'^) + xd{x, 



xu{x, fjr) + xd{x, iJ, ) 



(62) 



With the s and s distributions constrained to have the same number, jj^ slx^fi"^) dx — s{x, fi^) dx, Ag and A'^ are 
defined in terms of k, a and a' . In the CCFR analysis, four parameters, k, a, Aa = a — a' and the charm quark 
mass TOc, are determined from the data with th e charm hadron branching ratio fixed to the value obtained from other 
measurements, Bl = 0.099 ± 0.012 (see Section [ V E 6 ) . The results are: 



K = 0.536 ±0.030 ±0.037 



0.064 
0.098' 



a = -0.78 ± 0.40 ± 0.75 ± 0.98, 
Aa = -0.46 ± 0.42 ± 0.40 ± 0.65, 
rric = 1.66 ± 0.16 ± 0.07 t ^'^i GeV/c^ 



(63) 



where the first error is statistical, the second is systematic, and the third is due to the uncertainty in B].. The value 
of Aa = —0.46 ± 0.87 indicates that the momentum distributions of s and s are consistent, and the difference in the 
two distributions is limited to —1.9 < Aa < 1.0 at the 90% confidence level. Brodsky and Ma (1996) have compared 
their intrinsic strangeness model to these CCFR fits. The simple application of the model in their paper uses a scaling 
of the d valence quark distributions to do the QCD evolution of the intrinsic strange quark sea. The results of the 
calculation including a 30% extrinsic strange quark component are inconsistent with the CCFR analysis as shown in 
Fig. 4 of the Brodsky and Ma (1996) paper. From the figure, the predicted ratio is xs{x) / xs{x) = 1.7 at x = 0.07 
and falls to 1 at a; = 0.22. This difference corresponds to a Aa « —3 which is inconsistent with the CCFR limit. 

The CCFR fits assume a simple power law relation of the strange to non-strange sea and may not have sufficient 
flexibility to encompass an intrinsic strange quark contribution at high x. An additional high x component would show 
up in the CCFR analysis as a modification of B^ (and | VcdP) for the neutrino charm production and as an unexplained 
excess at high x for the antineutrino production. The agreement of the antineutrino data (Fig. |3^) and the agreement 
of the fit Be = 0.1091 ± 0.0097 with B] = 0.099 ± 0.012 (see Section |IVE6| ) extracted from other measurements 
limits the size of any intrinsic contribution to below 20%. Further analyses will be needed to directly test the intrinsic 
strangeness models more precisely. Since the intrinsic component mainly contributes at large x, it is unlikely that an 
s / s asymmetry could explain the discrepancy between the chargcd-lcpton and neutrino measurements of F2 in the 
small-x region. 



6. Measurements of \ Vcd \ 

If the CKM matrix elements are not assumed, then the CCFR NLO results (Bazarko et ai, 1995) can be interpreted 
in terms of a, mc and the following products : 

\Vcd\'^Bc = (5.34 ± 0.39 ± 0.24 t ^;25) x 10"^ (CCFR - NLO), 

|2 o /r, nn _i_nin _i_nnG+ 0.06\ . , 1 n— 2 



-iVas^B, = (2.00 ± 0.10 ± 0.06 1 ''^■'ll) X 10-^ (64) 

K + Z 

where the first error is statistical, the second systematic, and third from the QCD scale uncertainty. 



As shown in Table VI, there is little difference in the parameters if the analysis is performed to leading order or 
next-to-leading order. Accordingly, we may combine the CCFR result with the LO result from the CDHS collaboration 
(Abramowicz et al., 1982), 

= (4.1±0.7i°:J^) X 10"' (CDHS - LO), 

where the first error is the total experimental error and the second is the QCD scale error which is not given by the 
original analysis but assumed to be the same as in the CCFR measurement. Combining the two results, assuming 
that all of the experimental errors are uncorrelated, yields 

\Vcd\'^Bc = (4.92 ± 0.52) x 10"^ (CCFR/CDHS). (65) 
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These combinations can be used to extract |Vcdp and k|14sP when Be is determined from other data. Be is 
determined by combining the charmed particle semileptonic branching ratios measured at e~^e~ colhders (Barnett 
et al., 1996) with the neutrino-production fractions measured by the Fermilab E531 neutrino-emulsion experiment 
(Ushida et ai, 1988). Using an E^is > 30 GeV cut, E531 determined the following production fractions: 52 ±6% 
42 ± 6% 1 ± 2% and 5 ± 3% A+. In the E531 analysis, events that could not be unambiguously identified as 
_D+ or were all categorized as £>+ events. To remove this small bias, a re-analysis was performed that included 
updated values of the charmed hadron lifetimes (Bolton, 1994a; Bolton, 1994b). This re-analysis finds the following 
production fractions with an E^is > 30 GeV cut: 60 ± 6% 26 ± 6% D+, 7 ± 5% £>+, and 7 ± 4% A+. These 
production fractions are consistent with those measured by e+e~ experiments (Bortoletto et ai, 1988). With these 
particle fractions, B\ = 0.099 ± 0.012 (Bolton, 1994a; Bolton, 1994b) and when encorporated with the combined 
measurement of Eq. p| gives the value of the CKM matrix element 

\Vcd\ = 0.220 t mi (CCFR/CDHS), 

where the error indicates all sources of uncertainty including the fi^ scale uncertainty. This value compares well with 
the PDG value, \Vcd\ = 0.221 ± 0.003, which is determined from measurements of the other matrix elements and the 
unitarity constraint on the CKM matrix assuming three generations. The errors on the direct \Vcd\ measurement from 
V charm production are currently at the ±9% level which precludes a precise test of the CKM unitarity. 

The CKM parameter \Vcs\ requires an independent determination of the strange sea fraction, k, which is currently 
unavailable. Using the value of B], given above with the result of Eq. ^ and making the conservative assumption 
that K < 1 implies that 

\Vcs\ > 0.69 at 90% C.L. (CCFR). 

In the future, k could be determined from a measurement of the xF^ difference between neutrinos and antineutrinos, 
xF^{x, Q^) — xF^{x, Q^) — ix [s{x, Q^) — c{x, Q^)] ■ It should be possible with this technique to measure k to 20% 
which would allow the quantity |V^(i|/|Ks| to be determined to 10% from the dimuon measurements. 

F. Charm Sea Measurements from WSM Production 

WSM production can be used to probe the charm content in the nucleon through the neutral current scattering of 
a Vfj, off a charm quark which subsequently decays into a fj,~^ . A charm component is expected in QCD through gluon 
splitting into cc pairs. In a typical QCD calculation, the charm sea is set to zero for < Qg « 4 GeV^ and is then 
allowed to evolve via the DGLAP equations. In the interaction, NC scattering off a charm quark is suppressed due 
to the kinematic effects related to producing two heavy charm quarks (cc) in the final state. 

iy^ + c{+c) — > c (+c) (66) 

^ s + /i+ + i^^ 

For a mean energy of 100 GeV, this suppression introduces a factor of ~ 0.3 relative to massless quark production. In 
addition, to go from the charm production rate to the WSM rate, one must incorporate the semi-leptonic branching 
ration for the final state charmed particle (~ 10%) and the acceptance for the muon to be measured in the detector 
(typically ^ 20%). The size of the cross section for neutral current scattering off charm quarks in the nucleon can be 
estimated using the charm sea calculations. Using the NLO parton distributions of the CTEQ collaboration (Botts 
et at, 1993), the cross section for neutral current charm scattering relative to the total neutrino CC scattering cross 
section is a factor of 0.005 lower. Combining all these factors, one would expect a WSM rate which was approximately 
3 X 10~^ of the normal CC rate. 

The CCFR collaboration has made a study (Mishra et ai, 1989) of WSM production using the data from their 
narrow-band E616/E701 experiment. The dominant background comes from the contamination in the beam. 
CCFR minimizes this background by demanding the observed energy in the WSM events to have E^s > 100 GeV 
which reduces the background rate to 2.3 x 10~^. Other background sources (z/g induced dilepton production, mis- 
identified dimuon events, and NC interactions with a n/K decay in the hadron shower) contribute an additional 
background at the 1.5 x 10"'* level. For Eyis > 100 GeV, the CCFR experiment observes 43.0 ±6.6 WSM events with 
a calculated background of 26.9 ± 5.1 events leading to a 2a excess of 16.1 ± 8.3 events or a rate of (2.3 ± 1.2) x 10""* . 
The CCFR collaboration presents the excess as a 90% confidence level (CL) upper limit on the WSM rate of 4.3 x 10~^ 
which is much larger than the estimate of 1.4 x 10~^ for NC charm scattering with Eyig > 100 GeV. 
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V. ELECTROWEAK PHYSICS WITH NEUTRINO BEAMS 



Lepton scattering experiments permit the study of electroweak interactions at space-like momentum transfers of 
10~^ <Q^< 10+^ GeV^ in a fixed target environment with — > 10*^ GeV^ at the HERA ep collider possible in the 
future. As such they complement the large ^/s time-like measurements now performed with extraordinary precision 
at LEP, SLD, and the Tevatron. It is sometimes not appreciated that electroweak measurements performed in fixed 
target lepton scattering experiments can also still compete with the collider results in probing the Standard Model. 

Neutrino scattering, in particular, has contributed to our understanding of electroweak physics over roughly four 
historical phases since the early 1970's: (I) confirmation of the existence of neutral currents through first observation 
of y^e~^ Vfj_e~ at CERN (Hasert et ai, 1973); (II) early measurements of sin^ 9w that provided a critical ingredient 
to the successful prediction of the masses of the W and Z bosons (Arnison et ai, 1983); (III) second generation weak 
mixing angle determinations that probed one loop corrections to the Standard Model and provided among the first 
useful upper and lower limits on the top mass; and, finally, (IV) third generation experiments that seek to test the 
internal consistency of electroweak theory through precise coupling constant extrapolations that can be compared to 
other precision electroweak measurements at colliders. In this Review we will summarize the current experimental 
state of electroweak measurements with neutrino probes (period III) , and give short and intermediate term prospects 
for current and future experiments (period IV). Detailed discussions of results obtain prior to direct observation of 
the W and Z may be obtained elsewhere (Fogli, 1986; Kim et ai, 1981). The remainder of this section divides 
into four parts: Theoretical Motivation, Experimental Neutrino-Electron Scattering, Experimental Neutrino-Nucleon 
scattering, and Summary and Outlook. 



A. Theoretical Motivation 

Neutrino-nuclcon and neutrino-electron scattering depends on the left and right handed v — Z'^ and -target 
fermion couplings. Because all of the parameters of the Standard Model except the Higgs mass are now well known, 
the couplings measured in neutrino scattering are predicted to high accuracy. Any significant deviation between 
measurement and prediction would thus indicate new physics. Other experimental programs (LEP/SLD, LEP II, 
CDF/DO) also probe for new physics through precision electroweak tests, and neutrino experiments complement 
these high energy efforts along several lines: 

• By measuring different combinations of couplings than collider experiments, in particular, couplings to light 
quarks. 

• By measuring cross sections at moderate space- like momentum transfer, as opposed to the large time- like 
scattering explored at colliders. 

• By extracting cross sections with "orthogonal" analysis tools: for example, different radiative correction packages 
and very limited quark-gluon fragmentation model dependence. 

Electroweak measurements are frequently summarized by quoting a value of the weak mixing angle sin^ 9w however, 
this parameter can be defined in many different ways, e.g., 

1 - /Ml = sin^ Ow (on - shell) (67) 

and 



.^o^ _ (l - sin' Ow{cE))f - sinHw{e^) 



; \ - sin^ (eff ) ) + sin"* ew{e^) 



(68) 



with the differing definitions appropriate for different experiments. A less confusing way to compare the sensitivity 
of different experiments to electroweak physics is to translate all measurements to. an equivalent determination of a 
common physical parameter, for which a convenient choice is the W boson massQ. It follows simply that a ±0.002 



on-shell mixing angle measurement in neutrino scattering (see Section V A4 m) is equivalent to a ±100 MeV W^— mass 



measurement, which compares well with expected results from LEP II and CDF/DO. 



*These translations use SM electroweak theory with the fine structure constant a, the Fermi constant Gf, the Z° mass Mz, 
and the top quark mass Mtop, as input. Low energy electroweak observables can be calculated to sufficient accuracy without 
knowing the Higgs boson mass Mr- 
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1. Basic Cross Sections 



Muon neutrino scattering cross sections off light fcrmion targets can be expressed in tree-level as 



-2 



where / = e",it,d, s,c; Gp = (1.16639 ± 0.00002) x lO^^GcV"^ jg ^j^g Yeinii constant, Mz 91.187 ± 0.007 GeV 
is the mass; s is the effective center of mass energy, whic h dep ends on /; y is the inelasticity; and £f,rf are 



left and right handed coupling constants summarized in Table VIII. At fixed target energies, the propagator term 
(l + sy/M^) seldom differs from unity by more than one percent for quark targets and is entirely negligible for 
electron targets. Although all experiments apply propagator corrections, this factor will be omitted for brevity in 
subsequent formulas. 



./ 






e~ 
u, c 
d, s 


-\ -l-sin"' 6w 
i - 1 sin^ Ow 
-i + Isin^ Ow 


— 1 sin^ 9w 
i sin'^ 9w 



TABLE VIII. Neutral current standard model couplings for electron and quark targets. 
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Electron-neutrino-electron (t'eC" and v'eS~) scattering cross sections contain additional charged-current and inter- 
ference terms. These processes have been observed (Reines et ai, 1976; Allen et ai, 1990; Krakauer et ai, 1990; Allen 
et al., 1991; Krakauer et al., 1992a; Krakauer et ai, 1992b; Allen et ai, 1993); and experimental results have provided 
important qualitative confirmation of the SM, for example, the elimination of the four-fold sign ambiguity in axial 
and vector coupling constants of the electron. However, limited statistics preclude precision electroweak tests. Sim- 
ilar comments currently apply to very high NC scattering measurements at HERA (Derrick et al., 1995; Derrick 
et al., 1996b; Aid et al., 1996a); although there is considerable hope in this case for future improvement (See for 
example, Cashmore, Elsen, Kniehl, and Spicsbcrgcr, 1996). Accordingly, we will restrict our discussion to scattering 
experiments using i/^ and beams. 



2. Neutrino Electron Scattering 

For electrons, Sg = 2meEi, with rUe the electron mass and Ei, the neutrino energy. The total NC scattering cross 
sections then follow as 



dy 

da{v^e~ i^fie~) 
dy 



= cT-o'E^{el + rl{l-yf), (71) 
= a^,^E^{il{l-yr+rl), (72) 



with ctq^ = ^"^F™" = 17.2 X 10~^^ cm^/GeV. The electron mass sets the scale for the very small cross section and 
also drives the kinematics. With a maximum momentum transfer Q^^x = '^irieE^ ~ O.OOliJi/ GeV^, with E^ in GeV, 
the electron always scatters nearly parallel to the neutrino beam. 

The charged current process I'^e" — > t'e/U" , sometimes called inverse muon decay, represents an electroweak test 
in its own right. The tree level cross section is 

M.,e-^.eM-) ^ f ^ _ ^\ " 

dy 2meE,J ^ ' 

with the kinematic threshold occurring a.t Ei, = 11 GeV. The small center-of-mass energy forces the muon to be 
emitted in the forward direction. 

Muon decay measurements do not measiirc final state neutrino helicities and thus cannot distinguish between vector 
left-handed couplings, gY^, and scalar couplings ^f^; the inverse muon decay cross section is proportional to IqYl] ' 
and thus resolves this ambiguity. Muon decay and inverse muon decay are required to establish rigorously the V — A 
character of the charged current. 

Electroweak radiative corrections to i'^e~ and S'^e" NC (Novikov, Okun, and Vysotsky, 1993) and CC (Bardin 
and Dokuchacva, 1987) scattering have been calculated. An apparently coincidental cancellation between the two 
dominant radiative corrections in NC scattering causes the effective weak mixing angle sin^ 9'(y to be very close to 
sin^ 6w{eS) extracted from asymmetry measurements at LEP/SLD. 



3. Neutrino-Nucleon Scatteri,ng 

For quarks, the effective center-of-mass energy is Sq = 2ME,y^, with M the nucleon mass and ^ interpretable as 
the fraction of the nucleon's 4-momentum carried by struck quark. An experiment can only measure, in practice, 
the average over ^ of the sum over all quark targets inside the nucleon (and in fact, these distributions are further 
averaged over incident neutrino energy distribution): 

dai.,p-.u,X) ^2^^l^^^ E m + r){l-yf)^<l{^) 

^ q—u,d,s,c 

+ {ij(l-yf+rj)m)], (74) 
dai,,p^^,,X) ^2G%^J^^ E {[m-yr+r})m) 

^ q=u,d,s,c 

+ {ij+rjil-yf)^q{0]. (75) 
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Neutron cross sections follow fro m isos pin arguments with the substitutions u{u) — > d{d), d{d) u{u), s{s) s{s), 
and c(c) c(c) (but see Section V A 4 g). Contributions from b and t quarks are assumed to be negligible. 
For isoscalar targets, omitting second generation quark contributions for the moment, 

dajiy^N ^ i^^X) ^ GlME , ^^^^ 
dy TT 



m {rl + rl) (1 - yf) {^ui^ + Cd{0) 

+ {{€ + Q (1 - yf + rl + rl) (efi(e) + m))] , (76) 

^"^'^-^r ^-^^ - f [m + ^l) (1 - yf + rl + rl) i^uiO + m)) 

"2/ Jo 

+ {€ {rl + rl) (1 ~ y)') {^uiO + m))] ■ (77) 

At this level, one notes that NC cross sections can be re-expressed as 

Ty -{^u + e^} (78) 



,2 



r: 



2^^ d(j{9f,N n+X) 
dy 



da{v^N -> v^,X) _ I 2 , 2\ d(T{u^N ^ 11 X) 

dy -yru + r^) (ty) 



dy 

a result that follows more generally from isospin conservation (Llewellyn Smith, 1983). When integrated over y, the 
above imply the Llewellyn Smith (LS) relations 

= {el + il)+r{rl + rl), (81) 
i -sin2e'H' + (l + r)sin4 6'H'; (82) 

f ^ <ia(iJ^N-,v^X) 

Jdy '-'dy 

^{^l + Q+r-'{rl + rl), (84) 
^^-sin^0w + {l + r~^)sm'^ew; (85) 

with 

r = — ^ (86) 

Jdy '"dy 

the to Df^ CC charged current ratio. Because the LS relationships also hold at the differential cross section level, r 
can be defined to incorporate experimental hadron energy acceptance. For an ideal full acceptance experiment, r ~ 0.5, 
whereas typical experimental cuts reduce this value to r ~ 0.35 — 0.45. The LS formula is valid for an idealized isospin 
zero target composed only of first generation quarks. Corrections must be applied in real experiments to take into 
account the effects of non-isoscalar targets and, especially, heavy quarks involved in the scattering. 
The left handed coupling can be isolated through a linear combination of i?" and i?'': 

R- = (87) 



(88) 



d<j(v^,N~*v^X) _ dij{iJ^N~*v^X) 

dy dv 

da{u^,N^ti-X) _ drj{u^N^n+X) 
dy dy 

K + Q (89) 
^-sin^^w^. (90) 
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This final line expresses the Paschos-Wolfenstein(PW) relationship (Pashos and Wolfenstein, 1973). It retains its 
accuracy even if heavy quark contributions to the NC and CC cross section are included, provided that these con- 
tributions are the same for and i?^. This latter feature makes this quantity the most superior of all electroweak 
observables in the neutrino sector. 



4- SM Corrections to z/^(i/^)A'" Cross Sections 

Several corrections must be applied before SM couplings can be extracted from the LS or PW relations or their 
variants. 

a. SM Electroweak Radiative Corrections Radiative corrections approximately factor into two parts: QED final 
state radiation (De Rujula et al., 1979), for which final state muon bremsstrahlung diagrams are most important, and 
purely weak corrections which effectively shift coupling constants (Sirlin and Marciano, 1980). A complete treatment 
that combines the two effects exists (Bardin and Fedorenko, 1979). 

The effect of weak corrections can be represented schematically by the following: 

ef^{i + Ap)ef, (91) 

rf^{l + Ap)rf, (92) 
sin^ Ow {I + Ak) sin^ Ow (93) 

Calculation of Ap and Ak require a precise definition of s'rn^Ow- While any definition is possible, it has become 
conventional to use the on-shell or Sirlin definition (Eq. |6^) to describe lyN scattering. With this convention, the 
leading contributions to Ap and An which arise from the large top mass, 

largely cancel (Sirlin and Marciano, 1980) (Stuart, 1987); and a sin^ 6*1^ measurement from vN scattering can be 
thought of as effective VF— mass measurementO 

sin^ 9w{iyN) ~ sin^ 6\Y {on — shell). (96) 

Since the direct measurement of the top mass (Abe et ai, 1995b; Abachi et ai, 1995), the on-shell definition is less 
motivated, although it remains convenient so that the precision of neutrino measurements can be compared to that 
of colliders. 

QED radiative corrections are numerically large and must be handled with care. Their dominant effect is to harden 
the observed y distributions for Vf^N and Pfj^N CC cross sections relative to their expected theoretical forms; good 
control of acceptance corrections minimizes their effect. 

6. Structure Functions and Higher Order QCD Corrections The use of NC to CC cross section ratios coupled with 
the approximate validity of the LS and PW relations obviates the need to correct for light quark structure function 
effects. Similarly, uncertainties in as{Q^) contribute negligible error in the extraction of electroweak parameters. 

c. Heavy Quark Production The most uncertain parts of vN CC and NC cross section are due to scattering from 
and production of heavy flavor, 

I'f^s Vf_,s, (97) 
i^fiC -> Vf_,c, (98) 
iyf^d^p~c, (99) 
i^p.s — > /i^c, (100) 
I'/iS n^u; (101) 



^Strictly speaking, this is true only if sin^ 8w is extracted using the LS formalism, and even then, only if experimental 
acceptance effects do not move the effective value of r too far from r ~ 0.4 — 0.5. These restrictions have no consequence now 
that Mtop corrections can be made. 
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here s is counted as "heavy" . The phenomenology and experimental status of neutrino induced heavy quark production 
is described in Sec. [V. For the purpose of NC analyses, heavy flavor production introduces uncertainties in electroweak 
analyses through two main effects: 



• The neutral current processes are not simply related to analogous charged current cross sections through an LS 
type formula due to the role of the charm mass, CKM matrix elements, and large flavor asymmetry (s(^) ^ c(^)). 



• The threshold behavior of charm production is not known well enough, independent of the model used to 
describe it. In the naive leading order "slow rescaling" model where the threshold dependence in CC charm 
production is described by a single effective charm mass m^^^ , one only knows mf^^^ — 1.31 ± 0.24 GeV 
(Rabinowitz et ai, 1993), an 18% error. Better models of charm production exist, but these do not improve the 
experimental uncertainty. Threshold effects in charm production persist to very high neutrino energies because 
of the important role of sea quark mechanisms which predominate at low ^, where the effective center-of mass- 
energy 2AIEi,^ is not necessarily large compared to {nfj^^ . 

d. Isovector Effects Two kinds of isospin-violating effects are present; these can be thought of as nuclear and 
nucleon isospin violating effects. 

Heavy targets, such as iron, used in high statistics scattering experiments have a small (~ 6%) neutron excess. 
This excess creates cross section terms proportional to u — u + d — d = uy — dy and u — d, referred to as valence 
and sea contributions, respectively. The valence contribution has only directly been measured in relatively low 
statistics i^p and bubble chamber experiments (Jones et ai, 1994); however, global fits (see, for example, Lai et 
ai, 1996; Martin et ai, 1996) to parton distributions that compare high statistics i^N CC scattering to ep, fip, eD, 
and fj,D electroproduction provide tight indirect constraints. The sea quark contribution, now known not to vanish 
(Amadruz et al, 1992a), is not well-constrained, but is numerically small. 

Isospin violation exists even in an isoscalar target because electromagnetic and quark mass effects break the equality 
between u(proton) and d(neutron) quark distributions. These effects produce corrections to sin^ 9w of ^ 0.002 (Sather, 



e. Longitudinal Cross Section The simple decomposition of lyN cross sections into terms proportional to 1 and 
(1 — y)^ is broken by the existence of a non-zero longitudinal structure function Rl{^,Q^). As in the case of other 
higher order QCD effects, the longitudinal structure function modifies the y distribution but not the LS or PW 
relations. Influence on neutrino experiments arises only because the y distribution slightly affects CC/NC separa- 
tion experimentally. With sufficient care, these acceptance affects due to Rl{^,Q^) can be applied with negligible 
uncertainty added to coupling extractions. 

/. Higher Twist Effects Higher twist contributions to neutrino cross sections fall as l/Q^", with n > 1. This 
kinematic dependence alone severely suppresses their effects on NC analyses at high energies. Further suppression 
arises from the similarity of NC to CC cross section ratios from deep inelastic and higher twist sources, provided 
scattering occurs from an isoscalar target (a consequence of the LS relation). For example, purely elastic scattering, 
the "highest twist" of all, has a NC to CC cross section ratio of ~ 0.35 as compared to the DIS total of ~ 0.31. 

In fact, higher twist effects are arguably absent in any experiment that is in a regime where logarithmic QCD 
evolution of structure functions dominates (Q^ > 1 — 10 GeV^, depending on ^). In order to exhibit this approximate 
scaling behavior, the cross section must consist of a very large number of high twist terms; but, by duality arguments, 
many high twist terms added together are equivalent to the simpler quark scattering picture. 

Pumplin (1990) has argued to the contrary based on an application of vector/axial- vector dominance. In his model, 
the ai — p mass difference generates substantial and uncertain violations of the LS relationship (but not the PW 
relationship, since vector dominance produces equal i/^ and Dfj, cross sections). However, his model fails to describe 
CC differential cross sections unless its parameters are constrained to the level where effects on the NC to CC cross 
section ratios are small.Q Furthermore, from the arguments above, it may simply be incorrect to decompose an 



^Specifically, we have compared Pumplin's model to high precision charged lepton scattering and have concluded, conser- 
vatively, that his parameter So is constrained to So < 2 GeV^ at 90% confidence level. Pumplin's original analysis allowed 
So — > oo. 




(102) 



+2|v;.|'?s(0 + 2(i-y)'ec(0" • 



1992). 
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inclusive cross section displaying approximate scaling into a DIS contribution and a single higher twist term. More 
theoretical study of this issue would be useful. 

5. iVeui Physics Effects in i'p(!?^)e" and Ui_i{Ufj)N Cross Sections 

Direct effects of physics beyond the SM can appear in v^{v^)N scattering through higher order loop corrections (SM 
and non-SM Higgs, fourth generation quarks), new propagator effects (heavy , leptoquarks, compositeness), and 
lepton mixing (mirror fermions) (Langacker, Luo, and Mann, 1992). Somewhat surprisingly, the advantage enjoyed 
by neutrino measurements in probing certain classes of new physics derives often from the low energy scale of the 
process. For example, an unmixed new Z' would produce linear shifts in couplings in u^j^N NC scattering since both 
the Z^ and Z are far off-shell. In contrast, the effect of a heavy, unmixed Z' in e+e~ scattering on the Z^ pole would 
be to produce quadratic shifts in observed coupling since the Z' would not interfere directly with the Z^ at the pole. 
An advantage is also enjoyed over direct searches at the Tevatron, which have excellent sensitivity to constructive 
interference between and Z' , which increases lepton pair cross sections over expectations, but relatively poor 
sensitivity to destructive interference, which has the opposite effect. 

Other new physics not directly related to NC processes can manifest itself through shifts in the apparent values of 
R'^ , R'^ measured by experiment. For example, — *■ Vei^T oscillations would increase the value of R'^,R" relative 
to SM expectations because both and Vr charged current interaction produce an experimental NC signature 
(McFarland et ai, 1995). 

B. i^^e~ and i5^e~ Scattering: Experiment 

The experimental challenge in measuring NC processes is to separate genuine scattering from electrons from much 
larger cross section elastic and quasi-elastic and D^N reactions 

VeN ^ e~{n7r°)X, (103) 

where X consists of hadrons with too low energy to be measured; and from the rare, but kinematically similar, 
coherent pion production reactions: 

ly^N iy^TT°N. (104) 

Background rejection arises primarily from application of cuts that exploit the small momentum transfer inherent in 
scattering from electron targets. This demands good electron energy and angle measurement, which in turn mandates 
fine-grained low Z detectors. Since measurement of electron charge is not possible in a practical device, separate, 
high purity, sign-selected beams are required. Corrections must be applied for the wrong-sign neutrino content, and 
for the few percent component of v,, or in the beam. 

Charged current scattering is also optimally performed in a low Z, fine grained detector. This allows superior 
rejection of the primary background, y — > inclusive I'p^N scattering. At high energies, inverse muon decay can also 
be observed in dense detectors as the quasi-elastic processes become proportionally less important and the final state 
muon is easily observed. The low y v^N background processes can be measured from inclusive Pp,N scattering under 
the assumption that '^-i-^^r^t^' ^) = da(u,N^^+x) ^ y ^ 



1. Neutral Current Scattering 

CHARM II (Vilain et at, 1994a) recorded 2677±82 z^^e" and 2752 ±88 D^e~ events, after corrections, using a 690 
ton instrumented glass detector (Geiregat et ai, 1993) (Fig. ^ ) in a 2.1 x 10^^ proton-on-target exposure to the 450 
GeV CERN horn beam. The detector has excellent electron identification and measurement capabilities (Fig. ^ ). 
Average energies for i/^ and were 23.7 and 19.2 GeV, respectively. Backgrounds, mainly from coherent tt" production 
in elastic nuclear and nucleon scattering and quasi-elastic i^eN scattering, were at the 50% level. The CHARM II 
statistics ensure that this experiment dominates previous efforts at Brookhaven (Ahrcns et ai, 1983; Ahrens et al., 
1985; Abe et ai, 1986); hence, this section is essentially a resume of their results. 
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FIG. 43. The CHARM II detector. 
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FIG. 44. Display of two test beam events in CHARM II detector showing response to 10 GeV electron (a) and 10 GeV pion 
(b). The small dots represent tube hits and the area of the squares is a measure of energy deposition. 

CHARM II published their final results in terms of vector and axial vector couplings, gy' and g"^, which are simply 
related to left and right handed electron couplings, assuming the neutrino is purely left handed, by = ^ (ffy-^ + 97)' 
re = i (fif^^ - 5^^). The results are 

= -0.035 ± 0.012 ± 0.012 (CHARM II), 
g""/ = -0.503 ± 0.006 ± 0.016 (CHARM II), 

where the first uncertainty listed is the statistical error and the second is the systematic error. These measurements 
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were then combined to yield an effective weak mixing angle for ve scattering, sin , which, as noted earlier, is 
fortuitously close, theoretically, to the effective mixing angle measured at LEP and SLD: 

0.2324 ± 0.0058 ± 0.0059 (CHARM II), 



sm 



9ue 
W 



with the first error statistical and the second systematic. These results, based on an analysis of the absolute differential 
cross sections, were confirmed with an earlier study based only on the shape of the ^-distributions (Vilain et ai, 1993) 
(Fig. |45|). Excellent agreement with LEP/SLD results exists. 
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FIG. 45. Neutral current i^^ and j?,! differential cross section shapes measured by Charm II. 

The CHARM II data also provided, after combination with LEP data and lower energy ly^e^ scattering, an elegant 
demonstration of the flavor independence of the Z° coupling to neutrinos (Vilain et ai, 1994b). 

Finally, the comparison of CHARM II coupling measurements to SM electroweak theory allowed limits to be set on 
masses of new heavy Z' bosons. For the optimal case of no mixing with the ordinary Z'^ boson, the following limits 
hold (Vilain et ai, 1994c): 

Mz(^) > 262 GeV at 95% C.L. (CHARM II), 
Mzi^,) > 135 GeV at 95% C.L. (CHARM H), 
Mz(^) > 100 GeV at 95% C.L. (CHARM H), 
MziLR) > 253 GeV at 95% C.L. (CHARM H). 

More general limits are summarized in Fig. While these results have largely been superceded by direct searches 
and other precision measurements, it is nonetheless remarkable that such constraints follow from the modest statistics 
and low energy scale of CHARM II. 



64 




2. Inverse Muon Decay 

CCFR (FNAL E744/770) (Mishra et at, 1990) has a sample of 3248 ± 148 corrected events from an exposure 
corresponding to 3 x 10^ ordinary CC events at an average energy Ei, = 160 GeV. Charm II (Geiregat et ai, 
1990; Vilain et ai, 1996) recorded a signal of 15758 ±324 events from a much larger exposure (13.3 x 10^ CC events) 
at an average energy E^, — 23 GeV. Signal to background for this reaction was approximately 1 : 2.5 for both 
experiments. Figure shows the background-subtracted transverse momentum spectrum from CHARM II. 




0.02 0.04 0.06 0.08 0.1 

pf (GeV') 

FIG. 47. Distribution of inverse muon decay events as a function of squared transverse momentum pt from CHARM II. The 
solid line represents the expected distribution as simulated by a Monte Carlo calculation. 
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The — > oo asymptotic cross section slope (Eq. |7^) CTq'^ measured by the two experiments agree, and the total 
measurement errors are comparable. CCFR's error is dominated by statistics, whereas CHARM II, with their lower 
energy beam, is more affected by model uncertainties in making the background subtraction. The quoted cross section 
results, 

a^" = (16.93 ± 0.85 ± 0.52) x lO'^^cm^ (CCFR), 
a^" = (16.01 ± 0.33 ± 0.83) x lO^'^^cm^ (CHARM II), 

also agree with the SM prediction CTq'^ = 17.2 x 10^"^^ cm^. 

Consistency with the SM allows the two groups to set limits in the scalar couplings coupling of 

|(7fif < 0.30 at 90% C.L. (CCFR), 
l^fil^ < 0.475 at 90% C.L. (CHARM II); 

and CHARM II also quotes 

IgY^f > 0.881 at 90% C.L. (CHARM II). 

These results confirm the V — A character of the charged current interaction. 

It should be noted that CHARM II applied SM radiative corrections before setting their limits whereas CCFR 
apparently did not. The effect of the radiative correction was to increase CTq'^ by ~ 3%. Because the radiative 
corrections would (fortuitously) move the CCFR measurement closer to the SM, actual limits on scalar couplings are 
probably tighter. 



C. VfiN and i^^N Scattering: Experiment 



1. Method 



The critical issue in neutrino-nucleon NC to CC ratio measurements is the experimental separation between NC 
and CC events. This separation essentially depends on the presence or absence of a muon in a particular event, as 
can be seen in Figs. 48 and Ea, which show CC and NC candidate events in the CCFR detector. 



Run: 5153 Event: 1105 Date: Sat Dun 22 22:55:32 1996 
Triggers: nHHr*ir5ir6inDr9lNflllfl2lfl5l 




EHADO: 117.06 CeV 
EMUl: 26.40 GeV 




X-View 

FIG. 48. A charged current event candidate in the CCFR detector. 
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Two approaches have been successfully employed in high statistics measurements: event-by-event CC identification 
through muon track finding (CHARM (AUaby et ai, 1986; AUaby, et ai, 1987)) and statistical separation based on 
the distribution of longitudinal energy deposition in lAiN interactions [CDHS (Abramowicz et al., 1986; Abramowicz 
et ai, 1990), CCFR (Arroyo et ai, 1994)]. Figure ^ illustrates the latter idea as implemented by CCFR in the 
form of an event length distribution. The essential idea is that CC events produce penetrating muons that deposit 
energy over a large distance in the detector, whereas NC event lengths are characteristic of a hadron shower. CDHS 
employed a variant of this technique that allowed the separation length to vary with energy (Fig. ^). Corrections 
must be made for the fraction of short events which are actually CC, typically 5 — 20% and for the much smaller 
"punchthrough" of NC events into the long event length category. 
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FIG. 50. Event length distribution for CCFR neutral current analysis. Length is measured in units of scintillation counters 
exhibiting energy deposition in the event, with one counter corresponding to approximately 10cm of steel. The large peak at 
short event lengths consists mainly of neutral current interactions, with a background of Ue. NC and CC interactions and 
high y CC interactions. 



16000 



12000 



"i)! 8000 



4000 - 



NC 



CC 

Long . MONITOR 
Showers 

y 



1 1.4 



2 2.4 



FIG. 51. Scaled event length distribution for CDHS neutral current analysis. The shape of the plot is similar to that of 
CCFR, but the length separation is made energy depedent. 



The second most important experimental effect is the correction for v^, and in the beam. Both NC and CC 
electron neutrino interactions mimic the experimental signature of NC i/^ interactions, and their contribution to the 
apparent NC sample must therefore be subtracted. The subtraction is estimated by a detailed Monte Carlo simulation 
of the neutrino beam. This simulation can be tuned to describe l'e^ produced from charged kaon decay with high 
accuracy since the decay contribution is tightly constrained by measurements of the Vfi^v^ flux (Sec. ||). The 
largest uncertainty in the calculated electron neutrino flux comes from other sources of z/g such as beam scraping or 
decay. CCFR has recently directly measured the flux in their experiment using the difference in longitudinal 
energy deposition for electrons compared to hadrons (Romosan et al., 1996); the result confirms the calculated flux 
used in electroweak measurements. 

CHARM and CDHS attempted to extract values for and , and then infer sin^ 9w^ whereas CCFR varies the 
mixing angle in a garametric Monte Carlo approach to describe the data. Table IX compares properties of CCFR and 
CDHS/CHARM. El We focus on these three experiments in this section because other recent measurements (Bogert 
et al., 1985; Reutens et al., 1990) suffer from relatively poor statistical and systematic errors. 



Property 


CDHS/CHARM 


CCFR 


mean neutrino energy 


70 GeV 


160 GeV 


minimum hadron energy 


10 GeV 


30 GeV 


mean 


15 GeV^ 


35 GeV^ 


statistics 


2.0 X 10^ 


8.1 X 10^ 



TABLE IX. Comparison of properties of CCFR and CDHS experiments. 



''The latter pair ran in the same beam. 
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2. Standard Model Results 



The most accurate determination of sin ^6w from vN scattering is the recently updated final CCFR result (McFar- 
land et al, 1997): 

sin26'vK(on-shell) 0.2236 ± 0.0027cxpt. ± O.OOSOmodd (CCFR), 

where the first error is the total statistical and systematic experimental contribution and the second represents 
uncertainties due to model corrections. This value was extracted from the ratio of NC to CC events assuming values 
of m^f^ — 1.31 GeV, Mtop = 175 GeV, and Mh = 150 GeV; the explicit dependence on these parameters is 

sin26lvK(on-shell) = 0.2236 + 0.0111 (mf ^ - 1.31 GeV) 

+ (2.1 X 10-^) {Mtop - 175 GeV) - 0.0002 In (^-^^^^ (CCFR), 

where the masses are in units of GeV. The significant dependence of the result on the effective charm quark mass 
and the weak dependence on radiative corrections in translating to the on-shell mixing angle are evident. A detailed 
breakdown of systematic errors is provided in Table ^ The dominant contributions are from statistics, uncertainties 
in Ve production (mainly from neutral kaon sources), and charm quark production. 



SOURCE OF UNCERTAINTY 


S sin'^ 6w 


Data Statistics 


0.0019 


Monte Carlo Statistics 


0.0004 


TOTAL STATISTICS 


0.0019 


Charm Production (m^ = 1.31 ± 0.24GeV) 


0.0027 


Charm Sea 


0.0006 


Longitudinal Cross Section 


0.0008 


Higher Twist 


0.0010 


Non-Isoscalar Target 


0.0004 


Strange Sea 


0.0003 


Structure Functions 


0.0002 


Radiative Corrections 


0.0001 


TOTAL PHYSICS MODEL 


0.0030 


Ue Flux 


0.0015 


Transverse Vertex 


0.0004 


Energy Measurement 




Muon Energy Loss in Shower 


0.0003 


Muon Energy Scale (±1%) 


0.0004 


Hadron Energy Scale (±1%) 


0.0004 


Event Length 




Hadron Shower Length 


0.0007 


Vertex Determination 


0.0003 


Detector Noise, Efficiency 


0.0006 


Dimuon Production 


0.0003 


TOTAL EXPERIMENTAL SYSTEMATIC 


0.0027 


TOTAL UNCERTAINTY 




0.0041 



TABLE X. Uncertainties from updated CCFR extraction of sin^ 6w 
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CDHS and CHARM extracted their values for the weak mixing angle from the LS formula. CDHS (Abramowicz et 
at, 1990) obtains, using m^J^ = 1.50 GeV, Mtop = 60 GeV, and Mh = 100 GeV, 

sin26'M/(on-shen) = 0.228 ± O.OOScxpt ± 0.005,nodci (CDHS), 

and CHARM (AUaby, et ai, 1987) reports a very similar quality result 

sin 261^1/ (on-sheh) = 0.236 ± 0.005cxpt ± 0.005modoi (CHARM). 

Unfortunately these results were published at a time when the top mass was thought to be lighter; hence one must 
undo some of the radiative corrections to compare to the more recent CCFR result. If a common set of model 
parameters is employed, CCFR, CDHS, and CHARM agree in detail (Bolton, 1995), with the charm production 
uncertainty dominating in all cases. The average value of sin^ 9w obtained is 

sin^ ew{on shell) = 0.2256 ± 0.0035 (average vN). 

This value is obtained by fitting the CHARM, CDHS, and CCFR mixing angles and the CDHS and CCFR charm 
mass measurements to a simple model that allows sin^ 6w and m^^-^ to vary. The fit has a of 2.9 for 3 degrees of 
freedom and also yields an estimate of m^-''^ = 1.37 ± 0.19 GeV. The contribution to the sin^ 9w error from charm 
production is ±0.0023. 

Within the SM, the sin^ 9w measurement implies 

Mw = 80.25 ± 0.18 GeV (average i^N). 

The latter agrees well with direct determinations at the Tevatron and LEP II(The most recent high accuracy mea- 
surements are: Abachi et al., 1996; Abe et ai, 1995a; Ackerstaff et ai, 1996) Mw — 80.33 ±0.15 GeV, and the indirec^ 
determination from global fits to e+e" data at the Z^, (See, for example, Renton, 1995) Mw = 80.359 ±0.056 GeV. eI 

3. Model Independent Results 

The three collaborations have also attempted to present their results in a more model independent form, and we 
summarize these results here. Where possible, the explicit dependence of a parameter on the effective charm quark 
mass m^^f is given since correcting for CC charm production produces the largest systematic uncertainty in all results. 
In these cases, the charm mass contribution to the parameter is not included in the quoted error, but it can be easily 
calculated from m^J^ = 1.5 ± 0.3 GeV used by CHARM/CDHS and m^J^ = 1.31 ± 0.24 GeV used by CCFR. 

CHARM and CDHS quote values of R'^ and R" (Fig. [H), as well as separate extractions of Ap and sin^^vi^. 
CHARM obtains 

R" ^ 0.3093 ± 0.0031 (CHARM), 
R"" = 0.390 ± 0.014 (CHARM), 

using r = 0.456 ± 0.011 (CHARM); whereas CDHS extracts 

R'' = 0.3135 ± 0.0033 (CDHS), 
i?"" = 0.376 ± 0.0016 (CDHS), 

using r — 0.409 ± 0.014(CDHS). These cross section ratios represent the most model independent expression of i^N 
electroweak results; only corrections for experimental effects and the non-isoscalar nuclear targets are applied. As 
such, they cannot be directly used for electroweak tests without further model corrections. CDHS has applied these 
corrections to get 

R''^ = 0.3122 ± 0.0034expt - 0.009(m^^^ - 1.5 GeV) (CDHS), 
i?""" = 0.378 ± 0.014expt - 0.019(m^-'' ^ - 1.5 GeV) (CDHS). 



For a more detailed summary, see Barnett et al. (1996). 
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The "0" superscript denotes that these quantities have effects due to heavy quarks, radiative corrections, and other 
factors removed from the physical cross section ratios. Apphcation of these corrections requires shghtly different i/^ 

1.5 GeV) for neutrinos and 

expt ^ 0.004(m^// 



to CC ratios to be used in the LS formulas 



r" = 0.383 ± O.OMoxpt + 0.004(m^// 
1.5 GeV) for antineutrinos. 



0.42 



0.41 



0.40 



0.39 



0.38 



0.37 - 



0.36 




028 0.29 



0.33 



0.30 0.31 0.32 

FIG. 52. vs from CDHS. This measurement, along with a comparable result from CHARM imposed stringent con- 
straints on the top mass as soon as the Z boson mass was measured precisely. 



Because CDHS and CHARM took data in sign-selected beams, they could extract two electroweak parameters from 
their data without assuming the SM. Translated to the notation of Eq.'s 94-p5|, CDHS extracted values for these 
parameters (originally expressed as "p" and "sin^0vK")i using only their own data, of 

Ap = 0.009 ± 0.020cxpt - 0.023(m^^-^' - 1.5 GeV) (CDHS), 
Ksm^9w = 0.218 ±0.021oxpt - 0.011(m;;^-'' - 1.5 GeV) (CDHS). 

The larger experimental error on the mixing angle simply reflects the weaker constraints of a two-parameter fit. The 
strong top mass dependence of Ap and k in the SM allowed this measurement to be used to constrain Mt < 200 
GeV well before the direct observation of the top quark. CDHS and CHARM weak mixing angle measurements also 
decisively ruled out the once promising grand unified theory based on minimal SU(5) (which predicted sin^ 9\\r — 
0.218 ± 0.004), corroborating results from proton decay experiments (Sirlin and Marciano, 1980). 

CCFR quotes an experimental value (McFarland et ai, 1997) for a sum of couplings that is closest to what it 
actually measured in its mixed I'ti/i^ti beam: 

«= 1.7897 £2) 

-0.0916 {il - 1 

= 0.5820 ±0.0031 



-1.1479 (r^+r^ 
- 0.0782 {rl - 
0.0111(m;;^-^' - 



(105) 



1.31 GeV) (CCFR), 



This value can be used to constrain values of left- and right-handed couplings to quarks as shown in Fig. Using 
this measurement and the SM prediction of ksm = 0.5817 ± 0.0013, with SM parameters fixed from other precision 
electroweak measurements of Mz, Mw, and Mtop, CCFR sets limits on several new physics processes: 



71 



68% Confidence Allowed Regions for Couplings 




FIG. 53. Constraints on gj^ ~ iu+^d ^nd = r^ + r^ from CCFR model independent result and all neutrino data combined. 
Also shown are the projected constraints on these couplings from the NuTeV experiment. 

• For "left-left" four-fermion contact operator described by a Lagrangian — £ = ±{4tt / A^)lLj'^lLqLji,qL, 

> 4.7 TeV at 95% C.L., 
A7. > 5.1 TeV at 95% C.L. 



The approximate equality of these limits reflects the equal sensitivity to destructive and constructive interference 
effects in lyN scattering mentioned earlier and compares favorably with results set by CDF (Bodek, 1996). Limits 
on a set of other composite operators are also set. 

• For SU(5) Leptoquarks that do no induce flavor-changing neutral currents and which have dominant left-handed 
couplings, 

Mlq hiP^ > 0.8 TeV at 95%C.L. 



where M^q is the leptoquark mass and 77^ the left-handed coupling parameter. 

The CCFR result can be input into the formalism of Langacker, Luo, and Mawi (1992) to obtain parametric limits 
on a variety of new physics models. These results are summarized in Table Xl. El 



^The values in Table bC| are inferred from the CCFR measurement by the authors of this Review; the experiment may publish 
more extended results that differ slightly from those presented here. 
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Model 


95% 


C.L. Limit 


comment 




Mz, 


> 691 GeV 


Z{x) with 


it) =l^^'o^ 




Mz> 


> 215 GeV 


Z{x) with 


2|J =|sin='ew 




Mz> 


> 843 GeV 


Z{ip) with ( 


it) =§sin^ew 


Mz(i,),C = -\/f 


Mz- 


> 513 GeV 


/ 

Z{i>) with 1 


it) =§sin^ew 




Mz> 


> 54 GeV 


/ 

Z{Tp) with I 


it) =§sin^ew 


Mz(r,),C = 


Mz' 


> 101 GeV 


/ 

Z{r]) with ( 


=|sin2 0H^ 




Mz> 


> 879 GeV 


/ 

Z{r]) with ( 


it) =fsm^^w 




Mz> 


> 87 GeV 


Z{rj) with ( 


ft) =|sin^^w 


■Mz(3H),C = V s"* 


Mz- 


> 988 GeV 


Z{3R) with 


(it)^ = f-^^- 


M^fom, C = —A/la 


Mz- 


> 189 GeV 


Z{3R) with 


(|y = fsin^^^ 




Mz- 


> 244 GeV 


Z{3R) with 


(|)' = fsin^ew 


Non-SM Higgs 


Ap < 0.0049 


rad. corr. parameter 


SU(5) L(?pt()(|uai'k 


IMlq/vlI > 887 GcV 


RH coupling r/j? ~ 


extra i.iL foniiicjii 


siu- O'l < O.Dii 


mixing with 


(jixliuary iil 


extra fcririiori 


sin^ 


9^ < 0.037 


mixing with ordinary Uft 


extra cIl fermion 


sin^ 


9t < 0.009 


mixing with ordinary cIl 


extra, dn fermion 


sin^ 


9% < 0.074 


mixing with ordinary dn 


extra bl fermion 


sin^ 


91 < 0.043 


mixing with ordinary bl 


extra v^l fermion 


sin^ 


91' < 0.043 


mixing with ordinary f^L 


extra hl fermion 


sin^ 


9^ < 0.006 


mixing with ordinary fj,L 


extra i/^l fermion 


sin^ 


91" < 0.043 


mixing with ordinary Vf^L 


4-Fermi plus 


A+ > 4.1 TeV 


compositeness scale 


4-Fcniii minus 


A ; 


> 1.0 TcV 


compcjsitciK 


ss scale 



TABLE XI. Summary of limits on parameters for models of new physics beyond the SM as described by Langacker et al. 
based on CCFR preliminary model independent result. 
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VI. SUMMARY AND OUTLOOK 



High-energy, high-precision neutrino measurements have become an important tool for measuring the parameters 
of the Standard Model and for testing our understanding of the clcctrowcak and strong interactions. This Review has 
presented the current status of these measurements and has pointed out some of the open questions. These topics are 
summarized below along with comments on future measurements, now underway or being planned, which will have 
even higher precision and be able to address some of these questions. 

A. Structure Functions and QCD Tests 

Precision neutrino experiments have made important contributions to our knowledge of nucleon structure and the 
strong interaction through the measurements of the structtire functions and .TF3. Data have been taken over a 
wide range of targets and energies. The total neutrino cross section divided by energy has been found to be constant 
with energy, regardless of target. The structure function mcasin-emcnts have different assumptions concerning and 
radiative corrections, and are made on different types of target nuclei, making cross experiment comparisons difHcult. 
Iron data between CCFR and CDHSW can be directly compared and a discrepancy exists, particularly in the low-a; 
region and as a fimction of Q^. 

Neutrino scattering studies yield important determinations of the parton distribution functions. These results are 
important inputs to the global fits due to the unique flavor differentiation of the neutrino measurements combined 
with the good precision. The xF^i structure function is particularly important due to its close connection to the 
valence quark distribution. 

The recent neutrino measurements have comparable statistics to the charged-lepton structure function measure- 
ments. Comparisons of the neutrino F2 measurements to charged-lepton experiments allows an important test of the 
universality of the parton distributions measured in different processes. The comparison indicates general agreement 
over most of the kinematic region but a small discrepancy at low x after quark charge and nuclear corrections are 
applied. At present, this differences is not understood. The discrepancy may be related to assumptions concerning 
nuclear corrections or to systematic errors in one or more of the experiments at low x. The difference is too large to 
be explained by uncertainties in the strange sea which has been measured from neutrino charm production. 

Precise tests of QCD are now possible with neutrino scattering data. The CCFR E744/E770 data is the most 
precise and the evolution of the structure functions yields A^^~^^ = 337 ± 31(ea;p) MeV. This is equivalent to 
as{M'^) = 0.119 ± 0.002(ea;p) ± Q.004{scale). This is the most precise measurement of as at low and moderate Q"^ 
values. The agreement of the CCFR data with the QCD expectation is quite good. The GLS sum rule results have 
sufficient accuracy to test the predicted QCD corrections. Evaluation of this sum rule yields another measurement 
of the strong coupling constant which is consistent with the CCFR a,, result from QCD evolution. Because of the 
inclusive nature of the GLS and evolution measurements, these results have small theoretical errors and have been 
important in the past in establishing the QCD theory and more recently in making precise measurements of the 
fundamental parameters. 

The NuTeV experiment (Bolton et al., 1990) at Fermilab, which uses the upgraded CCFR detector and a new sign- 
selected beam (Bernstein et al, 1994), began taking data in May, 1996. This experiment will address many of the 
issues concerning structure functions raised in this review. Among the NuTeV goals are more precise measurements 
of as and R^. In order to investigate fully the discrepancy in the F2 measurements, more data in the low-x region 
with improved systematic understanding are required. The NuTeV experiment is expected to obtain approximately 
equal statistical errors on F2 with substantially reduced systematics. The most important systematic uncertainties 
are related to the absolute energy calibration of the detector. For the NuTeV experiment, precise hadron and muon 
energy calibrations will be obtained over a wide range of energies using precision test beam measurements throughout 
the run. With these improvements, the experiment should be able to measure the QCD scale parameter, Aqcd to 
better than 25 MeV and improve significantly the uncertainties in the F2 structure function measurements at low-a;. 

B. Measurements of the Strange Sea, Charm Sea, and \Vcd\ 

Neutrino charm production has been used to isolate the nucleon strange quark distributions, xs{x) and xs{x), and 
study the transition to the heavy charm quark. The CCFR NLO analysis agrees well with their dimuon data over a 
broad range of energies with a charm quark mass, rric = 1.70 ± 0.19 GeV. At NLO, the strange quark distribution is 
found to have the same shape as the non-stange sea but with a magnitude that is (48 ±5)% of a SU(3) symmetric sea. 
These measurements imply that an increased strange sea cannot be the explanation for the discrepancy between the 
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neutrino and charged-lcpton measurements of F2 at low-x. Investigations by the CCFR collaboration, also, indicate 
that the s{x) and 's{x) distributions are similar as expected from a perturbative QCD source but inconsistent with 
non-perturbative models of intrinsic strangeness in the nucleon. 

The CKM matrix elements \Vcd\ and \ Vcs\ are determined from the neutrino dimuon data using an effective charm 
semileptonic branching ratio found from i^-emulsion and e+e~ data. The combined CCFR/CDHS data yields a 
value of |Vcd| = 0.220 1 q'q2i which is in excellent agreement with the value of \Vcd\ = 0.221 ± 0.003 found from 
other measurements combined with the unitarity of the CKM matrix. With the assumption that the strange sea 
must be smaller than the non-strange sea (k < 1), |V^s| is limited by the dimuon data to be > 0.69 at 90% C.L.. 
The parameter k could possibly be determined in the future from a measurement of xF^{x,Q'^) — xF^{x,Q^) = 
4a; [s{x, Q"^) — c(x, Q^)] . A 20% measurement of k will provide a 10% measurement of • 

Present and future high-energy neutrino experiments should be able to significantly improve these measurements. 
The NuTcV experiment at Fermilab is running with a sign-selected beam which will eliminate the v/V confusion 
in the dimuon channel at low x with statistics similar to the CCFR (E744/770) sample. NOMAD (Astier et al., 
1991) and CHORUS (Armenise et al., 1990; DeJong et al, 1990) at CERN are presently running with a high-energy, 
horn-focused beam searching for neutrino oscillations. These experiments will also record a sizeable number of charm 
production events (about 20,000 dileptons for NOMAD and 20,000 reconstructed charm particles for CHORUS.) As 
in the E531 experiment, CHORUS will be able to reconstruct exclusive channels allowing precise tests of the theory. 
The ultimate neutrino charm production experiment is COSMOS(E803 at Fermilab) (Kodama et al., 1993). Like 
CHORUS, the experiment is a hybrid emulsion experiment that can reconstruct exclusive charm final states. The 
estimated data sample is 200,000 reconstructed charm events which might allow a determination of \Vcd\ to ^ 2% 
(Bolton, 1994a; Bolton, 1994b) . 

WSM production has also been used to probe the charm content in the nucleon through the neutral current 
scattering of a z/^ off a charm quark which subsequently decays into a /i"*" . From their WSM data with Eyis > 100 
GeV, the CCFR collaboration has set a 90% confidence level (CL) upper limit on the WSM rate of 4.3 x 10^'' which 
is much larger than the estimate of 1.4 x 10~^ for NC charm scattering. 

The currently running NuTeV experiment will have a substantially reduced contamination in their new SSQT 
beam which corresponds to a 2 x 10~^ WSM rate for Eyis > 100 GeV. If additional cuts can be used to minimize the 
other backgrounds, the NuTeV experiment may be able to isolate a NC charm scattering signal and make the first 
direct estimate of the charm sea in the nucleon. 



C. Electroweak Measurements with Neutrinos 

Experimental results on scattering on electrons and nucleons confirm predictions of the SM at low to moderate 
space-like momentum transfer, and in some cases, provides interesting limits on possible new physics beyond the SM. 

CHARM II, and earlier experiments, have demonstrated from NC measurements with electron targets that 
couplings to neutrinos are independent of family number and that the Z'^ couples to electrons in processes with very 
small space-like momentum transfer with precisely the strength specified by the SM and precision measurements at 
very high energy collider experiments. CCFR and CHARM II measurements of CC scattering of neutrinos from 
electrons demonstrate conclusively that the CC process is uniquely V-A and set tight constraints on possible exotic 
couplings beyond the SM. 

Precision NC measurements with nucleon targets by CDHS, CHARM, and CCFR allow an independent and com- 
petitive {AMw = ±190 MeV) prediction of the MK-boson mass within the SM in a very different kinematic regime 
than direct and indirect collider determinations. The consistency of neutrino measurements with SM predictions 
allows limits to be placed on a variety of possible new physics parameters, especially those involving light quarks. 

The only way to probe for the expected new physics beyond the SM is currently through precision electroweak tests, 
and this situation will likely persist at least until the LHC turns on and begins taking data. Because effects of new 
physics at low energies will be small and subtle, if they are there at all, it is desirable to perform as many cross checks 
as possible. Neutrino-nucleon scattering experiments are perhaps the only non-collider measurements that have the 
required sensitivity to provide these checks. And indeed, there are prospects for further gains from two near-term 
efforts that should improve electroweak measurements in the neutrino sector. 

The NuTeV experiment is currently taking data with high intensity sign-selected beams that will permit the first 
significant exploitation of the PW relationship, which will reduce the error due to heavy fiavor production by a factor 
of 4-5. The new beam design also eliminates much of the uncertainty attributed to if£ production of i/g- Depending 
in the total number of protons ultimately delivered, NuTeV should measure sin^ 6w with an error of 0.002 — 0.003, 
which is equivalent to a W^— mass measurement of 100-150 MeV. 
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High sensitivity searches for oscillations are now underway with CHORUS and NOMAD at CERN soon 

to be joined by COSMOS and MINOS (Abies et ai, 1995) at Fermilab. While these experiments probably will not 
directly measure electroweak parameters to interesting precision, they should greatly improve our understanding of 
neutrino induced charm production. It is possible that the error on the effective charm mass will be reduced by an 
order-of-magnitude or more, particularly by CHORUS and COSMOS, which will have unique abilities to measure 
inclusive charm production in their emulsion targets. With the largest systematic uncertainty in Vf^N perhaps nearly 
eliminated, one might envision a successor to the NuTeV experiment running at Fermilab with the higher intensity 
of the Fermilab Main Injector, or possibly even a new neutrino experiment at LHC or a muon collider. 
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